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SECTION  1 


INTRODUCTION 

The  work  discussed  here  is  an  experimental  means  to  estimate 
cloud  and  overcast  induced  probable  errors  in  a  geodetic  technique  of 
interest  to  AFTAC.  The  errors  of  interest  relate  to  the  difference  in  the 
assigned  times  for  the  initiation  of  a  surface  light  pulse  clocked  from  three 
wide  field-of-view  observing  stations  located  ir  space.  The  errors  are 
related  to  the  degree  of  time  smearing  a  fast  light  pulse  suffers  due  to  the 
multitude  of  path  lengths  taken  by  the  photons  which  contribute  to  the  re¬ 
ceived  pulse.  This  report  presents  experimentally  derived  estimates  of 
the  pulse  shape  one  would  expect  from  a  delta  function  light  source  at  or 
near  the  earths  surface  at  a  wide  field-of-view  optical  receiver  in  earth 
orbit  after  transit  through  a  sampling  of  natural  clouds.  The  experimental 
technique  employed  does  not  attempt  a  direct  measurement  of  the  shape 
of  the  received  light  pulse  (a  more  valid  but  more  costly  method)  but  rather 
infers  the  pulse  shape  from  cloud  geometric  and  radiance  measurements 
and  radiation  transport  theory. 

The  work  covered  by  this  report  was  supported  by  the  Air  Force 
Technical  Applications  Center  (TFE),  Patrick  AFB,  FL.  The  experimental 
phase  of  the  work  was  conducted  at  Kennedy  Space  Center  during  the  summer 
of  1977  in  conjunction  with  the  Thunderstorm  Research  International  Program, 
Trip  77.  The  data  reduction  and  analysis  phase  has  been  carried  out  at 
HSS  Inc,  Bedford,  MA  using  the  computer  located  at  the  Air  Force  Geophysics 
Laboratory,  Hanscom  AFB.  This  is  the  final  report  for  this  project  and 
includes  revised  versions  of  material  presented  earlier  by  interim  report 
and  by  papers  read  at  the  Satellite  Working  Group  Meetings  held  at  Patrick  AFP. 
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SECTION  2 


EXPERIMENTAL  CONCEPT 

2.  1  ELEMENTS  OF  THE  PROBLEM 

The  presence  of  an  overcast  or  clouds  in  the  line-of-sight 
degrades  the  capability  of  the  geodetic  technique  based  on  the  time- 
of-arrival  of  a  sharply  rising  light  pulse  at  three  or  more  stations  with 
wide  field-of-view  receivers.  The  loss  in  capability  arises  directly 
from  the  reduction  in  the  rate-of-rise  of  the  leading  edge  of  the  light 
pulse  with  a  consequent  loss  in  resolution  for  the  time  assigned  to  the 
arrival  of  the  pulse.  The  degree  to  which  such  a  system  would  be  de¬ 
graded  can  be  estimated  from  the  pulse  shape  one  would  expect  from  a 
delta  function  light  pulse  after  transit  through  various  cloud  types.  The 
experimental  concept  presented  here  was  developed  to  provide  a  means  of 
extracting  such  pulse  shapes  from  experimental  measurements  on  natural 
clouds  and  to  provide  data  useful  in  estimating  the  geodetic  errors  which 
would  prevail  for  the  cloud  types  measured. 

Any  photon  originating  in  the  delta  function  light  pulse  which 
transits  from  source  to  receiver  along  the  direct  path  arrives  at  the  re¬ 
ceiver  in  the  minimum  elapsed  time  after  the  occurrence  of  the  pulse.  A 
photon  traveling  by  any  other  path  suffers  a  delay  compared  to  the  minimum 
time,  the  magnitude  of  the  delay  depending  on  the  extra  path  length.  In 
clear  weather  essentially  all  of  the  energy  arriving  at  the  receiver  travels 
by  the  direct  route.  In  the  case  of  an  overcast  or  intervening  clouds  only 
a  tiny  fraction  of  the  energy  arrives  by  the  direct  route  and  most  arrives 
by  a  multitude  of  paths  dictated  by  the  scattering  conditions  in  the  intervening 
cloud. 


Figure  2. 1  presents  three  possible  photon  paths  between  an 
isotropic  source  near  the  ground  and  a  receiver  located  in  space  (the 
direct  path  from  source  to  receiver  is  not  included).  Paths  of  the  type 
shown  by  A  result  from  a  single  photon  scatter  by  an  air  molecule  or 
aerosol  particle,  the  scattering  particle  may  or  may  not  be  located  in  an 
identifiable  cloud.  If  the  location  of  the  scattering  particle  can  be  estab¬ 
lished,  the  total  path  length  can  be  computed  and  the  excess  over  the  direct 
path  determined.  This  is  a  purely  geometric  problem  and  leads  to  geo¬ 
metric  delay.  Paths  of  the  type  shown  by  B  result  from  multiple  scattering 
of  photons  along  the  direct  path  where  the  emergent  photon  continues  to 
travel  along  the  direct  path  from  source  to  receiver.  The  extra  path  length 
leads  to  scattering  delay  and  results  entirely  from  the  multiple  scatter 
with  the  geometry  of  the  source  and  receiver  playing  no  part.  A  more  probable 
path  is  shown  by  C  in  which  the  extra  path  length  results  from  both  geometry 
and  scattering  and  leads  to  a  combined  geometric  and  scattering  delay.  The 
received  pulse  shape  is  the  sum  over  time  of  all  photons  arriving  with  their 
respective  geometric  and  scattering  delays. 

The  most  direct  solution  to  the  problem  of  acquiring  cloud  mod¬ 
ified  pulse  shapes  is  one  of  direct  measurement  with  a  suitable  source  and 
receiver  deployed  as  shown  in  Figure  2.  1.  In  practice  this  approach  is  both 
difficult  and  costly.  An  indirect  solution  to  the  problem  is  presented  here 
which  employs  pulse  shape  synthesis  from  three  quantities.  For  every  line- 
of-sight  contributing  to  the  final  pulse  shape  these  quantities  are: 

1.  The  probability  that  a  photon  will  arrive 
along  that  path 

2.  The  geometric  delay  for  that  path 

3.  The  mean  scattering  delay  for  that  path 
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The  experimental  concept  now  becomes  one  of  devising  a  means 
of  recording  cloud  data  from  which  these  quantities  can  be  derived.  The 
remainder  of  this  section  is  devoted  to  a  review  of  experimental  concept 
and  a  non-mathematical  presentation  of  how  the  expected  pulse  shape  is 
synthesized.  A  full  development  of  the  system  is  given  in  subsequent  sections. 

2.2  EXPERIMENTAL  APPROACH 

The  geodetic  measurement  application  of  primary  importance 
concerns  pulsed  isotropic  sources  at  or  near  the  surface  with  the  receivers 
in  space.  The  approach  presented  here  make  two  fundamental  transform¬ 
ations  in  the  problem  for  experimental  reasons:  (1)  the  receiver  is  placed 
on  the  ground,  and  (2)  a  steady,  as  opposed  to  a  pulsed,  light  source  is 
employed,  the  sun.  These  transformations  permit  a  valid  and  experimentally 
simple  approach  to  the  problem.  The  first  transformation  is  based  on  the 
reciprocity  theorem  (discussed  fully  below)  which  states  that,  when  only 
linear  processes  are  involved,  the  source  and  receiver  can  be  interchanged 
with  no  change  in  result.  The  second  transformation  is  based  on  the  ability 
to  extract  from  the  s  teady  state  experimental  data  the  geometric  and 
scattering  delays  required  to  synthesize  the  pulse  shape. 

The  primary  experimental  procedure  consisted  of  photographically 
recording  with  all-sky  cameras  the  radiance  of  the  bottom  of  sun  illuminated 
clouds  simultaneously  in  four  wavelength  bands:  red,  blue,  narrow  band  at 

o  o  o 

7606  A,  and  narrow  band  at  7530  A.  The  narrow  band  at  7606  A  is  located 
in  one  of  the  Og  absorption  bands  of  the  atmosphere  (ON-BAND)  and  the 
narrow  band  at  7530  A  is  adjacent  but  not  in  the  absorption  band  (OFF-BAND). 
Each  of  these  photographs  provide  data  for  all  azimuth  angles  around  the 
recording  station  and  for  all  zenith  angles  down  to  the  horizon.  The  exposed 
photographic  film  was  processed  under  sensitometric  control  to  permit  the 
extraction  of  cloud  bottom  radiance  in  absolute  radiometric  terms.  In 


addition  to  the  primary  cameras  stereo  photographs  of  the  clouds  were 
recorded  from  two  separated  stations  by  cameras  with  wide  angle  lenses. 
Various  calibration  and  support  equipment  completed  the  experimental 
plan.  The  applicability  of  these  choices  and  the  manner  by  which  this  plan 
supplies  the  data  required  to  synthesize  an  expected  pulse  shape  can  be 
summarized  by  the  following  statements. 

1.  The  sun  has  a  spectral  power  distribution  similar 
to  that  of  the  sources  of  interest  and  is  therefore 
an  acceptable  substitute  for  these  sources. 

2.  Digitization  of  the  four  primary  photographs  pro¬ 
vides  data  for  all  lines -of-sight  of  interest  where 
each  data  element  corresponds  to  one  line-of-sight. 

3.  Registration  of  the  four  primary  photographs  brings 
all  data  elements  of  the  photograph  from  each  wave¬ 
length  band  into  coincidence  which  allows  any  one 
data  element  from  all  four  photographs  to  be  assigned 
to  a  single  line-of-sight. 

4.  The  stereo  pair  photographs  provide  the  height  of  the 
cloud  bottoms. 

5.  The  absolute  radiance  along  any  line-of-sight  is  a 
measure  of  the  transmission  from  sun  to  camera  and 
can  be  related  to  the  probability  that  a  photon  will 
arrive  along  that  line-of-sight. 
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6.  Any  line-of-sight  which  is  associated  with  an  area 
of  blue  sky  and  does  not  terminate  on  the  bottom 
of  an  identifiable  cloud  can  be  located  from  the 
ratio  of  the  blue  radiance  to  the  red  radiance.  Such 
lines -of-sight  contribute  a  negligible  fraction  of  the 
total  energy  received  and  will  be  ignored  during  data 
processing. 

7.  A  portion  of  the  total  geometric  delay  can  be  deter¬ 
mined  from  the  direction  of  any  line-of-sight  and 
the  cloud  bottom  height.  This  will  be  discussed 
more  fully  below. 

8.  The  scattering  delay  along  a  line-of-sight  due  to 
photon  diffusion  within  the  cloud  can  be  evaluated  by 
comparing  the  radiance  in  the  narrow  band  at  7  606  A 

o 

(ON-BAND)  to  radiance  in  the  narrow  band  at  7530  A 
(OFF-BAND).  Any  ON-BAND  radiance  value  is  re¬ 
lated  to  the  total  path  through  the  atmosphere  that  is, 

(1)  from  the  top  of  the  atmosphere  to  the  cloud  top, 

(2)  through  the  cloud,  and  (3)  from  cloud  bottom  to  the 
camera.  This  will  be  discussed  more  fully  below. 

The  equipment  to  implement  this  experimental  plan  was  designed, 
constructed,  and  operated  at  Kennedy  Space  Center  during  the  summer  of 
1977.  A  more  complete  description  of  the  equipment  is  given  in  Appendix  A. 

2.  3  THE  RECIPROCITY  THEOREM 

The  reciprocity  theorem  concerns  the  reversibility  of  photon 
propagation.  The  theorem  was  first  stated  by  Helmholtz  (Ref.  1)  and,  from  an 
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experimentalist's  point  of  view,  the  theorem  is  used  in  a  most  articulate 
elaboration  by  Minnaert  ,  (Ref  2).  Helmholtz's  original  statement  is 
translated  here  as: 

"From  point  A  a  unit  quantity  of  light  of  a  given  color, 
polarized  in  a  given  direction  a,  goes  in  such  a  direction  that 
after  a  series  of  reflections  and  refractions  the  quantity  x  finally 
arrives  at  B  and  is  polarized  in  the  direction  /3  .  If  we  now  allow 
a  unit  quantity  of  j8  polarized  light  to  start  out  backwards  in  the 
reverse  direction  from  this  final  ray,  then  after  the  light  under¬ 
goes  all  of  the  former  reciprocal  and  reversible  processes  the 
same  quantity  x  and  the  same  fraction  of  the  a  polarized  light 
arrives  at  the  initial  point.  " 

Minnaerts1  comments  are  significant.  "This  principle  is  widely  applicable: 
it  applies  to  reflecting,  refracting,  absorbing  and  scattering  media,  for 
each  separate  wavelength  and  for  each  plane  of  polarization.  It  fails  only 
in  the  cases  of  fluorescence  and  of  magnetic  rotation.  "  If  Minnaert  were 
writing  now  he  might  have  described  the  principle  as  one  including  all 
linear  processes  and  excluding  only  non-linear  processes. 

The  reciprocity  theorem  was  written  to  describe  source  and 
receiver  exchanges  in  steady  state  situations.  To  apply  it  to  the  problem 
addressed  here  requires  that  the  theorem  be  expanded  to  cover  transient 
light  signals.  The  authors  of  this  report  are  unaware  of  any  publications 
which  discuss  a  time  dependent  reciprocity  theorem  although  some  techniques 
of  Monte  Carlo  calculations  imply  time  dependent  reciprocity.  Two  state¬ 
ments  can  be  developed  from  the  reciprocity  theorem  (1  and  2  below)  and  a 
"reasonable"  additional  statement  is  added  (3  below)  which  implies  a  time 
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dependent  reciprocity  theorem.  The  validity  of  time  dependent  reciprocity 
can  be  demonstrated  for  specific  cases  but  no  general  theorem  has  been 
found  or  developed  as  part  of  the  work  reported  here.  The  three  statements 
are: 

1.  In  the  presence  of  any  condition  of  weather  or  of  ground 
reflectivity  the  power  received  from  the  sun  by  an  isotropic  detector  is 
in  fixed  proportion  to  the  power  which  would  be  received  at  the  same 
time  by  a  satellite  in  the  direction  of  the  sun  from  an  isotropic  source  at 
the  same  location  as  the  isotropic  detector.  This  assumes  that  the  spectral 
power  distribution  of  the  isotropic  source  corresponds  to  that  of  the  sun. 

2.  The  isotropic  detector  of  statement  1  above  can  be  re¬ 
placed  in  the  argument  by  a  receiver  having  a  specific  solid  angle  of 
response  if  only  the  isotropic  source  is  replaced  by  a  source  radiating 
into  the  same  solid  angle  and  direction. 

3.  Relative  to  t  ,  the  direct  passage  transit  time  of  a  photon 
from  sun  to  the  location  of  the  receiver  of  statements  1  or  2  above,  the 
probability  that  a  photon  will  arrive  at  some  later  time  t,  |3  (t-tQ)  is  the 
same  as  the  probability  of  a  photon's  delayed  arrival  at  the  satellite 
/3(t-tj)  when  the  appropriate  source  is  substituted  for  the  receiver.  Here 
tj  jLs  the  transit  time  for  direct  passage  from  the  location  of  the  ground 
based  source  to  the  satellite. 

A  significant  conclusion  is  reached  from  statement  1  above. 

That  is:  for  any  weather  conditions  during  the  day,  the  radiant  power  from 
the  sun  measured  by  an  isotropic  detector  in  any  location  can  be  used  to 
calculate  the  transmission  of  power  from  an  isotropic  source  at  that  loca¬ 
tion  to  a  satellite  in  the  direction  of  the  sun. 
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Statement  3  above  indicates  in  part  that  if  the  probability  of 
time  delay  of  photons  traveling  from  the  sun  to  an  isotropic  detector  is 
known  then  the  probability  of  delay  from  an  isotropic  source  to  a  detector 
in  the  direction  of  the  sun  is  also  known.  It  follows  that  a  pulse  shape 
derived  from  a  ground  based  detector  can  be  used  to  infer  the  pulse  shape 
which  would  have  been  measured  at  a  satellite. 

2.4  GEOMETRIC  DELAY 

Referring  again  to  Figure  2. 1  photon  paths  of  the  type  illustrated 
by  A  are  not  a  matter  of  concern  in  the  study  reported  here.  Such  paths 
occur  in  areas  outside  an  identifiable  cloud;  that  is,  in  blue  sky  areas,  and 
are  considered  negligible  contributors  to  the  total  received  energy.  Paths 
of  the  type  illustrated  by  B  are  special  cases  leaving  paths  of  the  type 
illustrated  by  C  for  consideration. 

The  photon  path  illustrated  by  C  in  Figure  2.  1  is  one  of  a  family 
of  possible  paths  from  the  bottom  to  the  top  of  the  cloud  (assuming  the 
source  on  the  ground  and  the  receiver  in  space).  The  one  shown  leaves 
the  top  of  the  cloud  at  a  point  directly  above  and  a  distance  B  from  the 
point  of  entry  at  the  bottom  of  the  cloud.  This  path  is  representative  of 
the  average  geometric  delay.  The  extra  path  length  which  leads  to  delay 
is  the  difference  between  the  direct  path  from  source  to  receiver  and  the 
path  consisting  of  (1)  source  to  cloud  bottom,  (2)  the  distance,  B,  and  (3) 
top  of  the  cloud  to  receiver.  The  data  reduction  procedures  used  this 
definition  of  geometric  delay. 

In  the  absence  of  absorption  the  transmission  of  a  cloud  or 
overcast  is  a  function  of  the  ratio  of  the  cloud  thickness,  B,  to  the 
scattering  mean  free  path,  jt  .  That  is  the  absorption  free  transmission 
is  a  function  of  B Ijl  .  This  ratio  can  be  calculated  for  each  line  of  sight 
by  proper  data  reduction  using  the  absolute  radiance  of  the  bottom  of  the 


cloud  from  the  red  photograph  taken  in  the  wavelength  region  X  =  6700  A 
(see  subsequent  sections). 

This  method  of  deriving  the  cloud  thickness,  B,  from  the 
recorded  data  makes  the  experimental  plan  given  above  complete  with 
no  independent  measurement  of  cloud  thickness  (or  cloud  top  heights) 
required.  Once  data  reduction  was  underway,  this  method  proved  inadequate 
for  a  significant  fraction  of  the  lines -of-sight  leading  to  spurious  results 
for  these  data  elements.  The  computer  programming  automatically  ignored 
such  results  and  excluded  the  data  point  from  contributing  to  the  final 
pulse  shape.  In  some  severe  cases  a  third  of  the  total  data  was  being 
rejected  by  this  criterion.  At  this  point,  a  modification  was  made  in 
the  data  procedures  to  permit  the  quantity  B  to  be  assigned  on  the  basis  of 
typical  published  values  for  the  cloud  types  being  considered.  With  this 
change  only  a  few  per  cent  of  the  lines -of-sight  were  automatically  rejected. 

2.5  SCATTERING  DELAY 

The  photon  path  marked  C  in  Figure  2.  1  includes  scattering 
delay  in  addition  to  geometric  delay  due  to  the  random  path  followed  by  the 
photon  in  its  transit  from  the  bottom  to  the  top  of  the  cloud.  The  random 
path  increases  the  total  path  and  the  increase  over  the  distance  B,  the 
cloud  thickness,  leads  to  the  scattering  delay.  The  extra  path  length,  and 
therefore  the  scattering  delays  which  occur  due  to  light  diffusion  through 
a  cloud  may  be  evaluated  by  correlating  the  strength  of  an  observable 
optical  absorption  structure  of  the  atmosphere  to  the  path  length  required 
to  produce  such  absorption  strength.  Part  of  this  path  length  is  due  to 
light  paths  external  to  the  cloud  which  can  be  calculated  as  in  the  case  of 
the  geometric  dealy.  The  remainder  of  the  absorption  related  path  length 
is  then  due  to  the  average  total  random  walk  path  length  of  the  photons 
traversing  the  cloud.  The  average  excess  path  length  is  related  to  the 
average  scattering  delay  through  the  velocity  of  light.  For  this  study  the 
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atmospheric  absorption  feature  measured  was  that  of  molecular 

o 

oxygen  at  7606  A  ,  the  ON-BAND  camera. 

As  stated  in  the  previous  section  in  the  absence  of  absorption 

the  transmission  of  a  cloud  or  overcast  is  a  function  of  the  ratio  of  cloud 

thickness,  B,  to  scattering  mean  free  path  /;  that  is,  to  B Ijt  .  When 

absorption  is  present  the  ratio  of  cloud  transmission  with  absorption  to 

2 

transmission  without  absorption  is  a  function  of  B  // .  This  latter  ratio  can 
be  calculated  for  each  line-of-sight  to  the  cloud  by  the  proper  radiometric 

C 

reduction  of  the  pair  of  photographs  taken  at  wavelengths  X  =  7606  A 

0 

and  X  =  7530  A.  the  ON-BAND  and  OFF-BAND  pair  of  photographs. 

2.  6  PULSE  SHAPE  SYNTHESIS 

The  procedures  and  quantities  outlined  above  were  used  to 
derive  the  pulse  shape  of  power  versus  time  one  would  expect  at  a  detector 
in  space  from  an  isotropic  delta  function  light  source  near  the  earth's  surface 
through  the  cloud  conditions  which  prevailed  at  the  time  the  experimental 
measurements  were  made.  These  predicted  pulse  shapes  are  the  final 
output  from  this  study.  The  synthesis  was  accomplished  by  computer  pro¬ 
cessing  the  experimental  data  in  the  following  manner. 

Assuming  all  calibration  and  other  preliminary  data  processing 
steps  completed,  each  pulse  shape  synthesis  starts  with  the  following 
quantities. 

1.  Solar  irradiance  on  the  top  of  the  cloud  or  overcast. 

2.  Azimuth  and  zenith  angle  to  the  sun. 

3.  Height  of  the  cloud  bottom. 

4.  Height  of  the  cloud  top  (except  when  this  quantity  is 
to  be  derived  from  the  data). 

5.  The  azimuth  and  zenith  angle  for  every  line-of-sight 
which  contributes  to  the  received  pulse  (this  excludes 
all  lines -of-sight  that  terminate  in  blue  sky). 
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6.  The  absolute  radiance  of  the  cloud  bottom  in 
the  red  band  along  every  contributing  line -of  - 
sight. 

7.  The  ratio  of  the  radiance  of  the  cloud  bottom  for 
every  contributing  line-of-sight  in  the  oxygen 
absorption  band  (ON-BAND)  and  just  outside  the 
oxygen  absorption  band  (OFF-BAND). 

By  the  methods  outlined  above  and  explained  in  detail  in 
subsequent  sections  the  following  quantities  are  derived  for  each  line-of- 
sight  that  contributes  to  the  final  pulse  shape. 

1.  The  rectangular  coordinates  of  the  tops  of  the 
clouds  corresponding  to  each  line-of-sight  to  the 
bottom  of  the  cloud. 

2.  The  mean  transit  time  from  the  top  of  the  cloud  to 
the  point  of  detection  along  each  line-of-sight.  (This 
statement  refers  to  the  measurement  configuration 
with  sun  and  ground  based  detector. ) 

3.  The  transmission  of  the  cloud  along  each  line  of 
sight. 

The  reciprocity  theorem  is  applied  to  the  interpretation  of  the 
values  of  transmission  and  mean  time  of  transit.  In  the  process  outlined 
so  far  values  of  transmission  and  transit  time  have  been  calculated  in  terms 
of  downward  flowing  light  from  the  sun  (the  source)  to  the  isotropic  receiver 
(the  camera  systems).  The  reciprocity  theorem  allows  these  values  to 
be  used  in  treating  a  delta  function  isotropic  source  located  at  the  position 
of  the  measuring  equipment  and  a  receiver  in  the  direction  of  the  sun. 
Corresponding  to  each  of  the  equal  elements  of  solid  angle  of  the  camera 
system  there  are  equal  solid  angle  elements  of  the  source.  When  pulsed, 
equal  amounts  of  energy  are  invested  in  each  element  of  solid  angle. 


Reciprocity  theorem  allows  assignment  for  each  of  these  elements  of 
solid  angle  a  mean  transit  time  and  a  transmission  value  to  the  corres¬ 
ponding  element  of  area  at  the  top  of  the  cloud. 

Let  the  number  of  data  elements  of  the  all  sky  camera  picture 
be  N  with  each  element  having  an  unique  assigned  number  n.  (1  <  n  <  N). 

If  the  transmission  associated  with  element  n  is  T  then  for  those  elements 

n 

exposed  to  blue  sky,  0.  Let  the  rectangular  coordinates  of  the  location 

of  the  cameras  be  x  -  y  -  z  =  0  then  for  each  of  the  elements  exposed  to 

clouds  there  is  an  element  of  area  at  the  top  of  a  cloud  with  coordinates 

x  ,  y  ,  z  ,  with  associated  transmission  T  and  associated  mean  transit 
n  •'n  n  n 

time  from  (0,0,0)  of  t  .  If  a  satellite  is  located  in  the  direction  of  the 

n 

sun  at  coordinates  x  ,  y  ,  z  then  the  time  of  flight  of  light  from  the  source 

s  s  s 

directly  to  the  satellite  is  t  where 


V 


2  2  2 
x  +  y  +  z 
s  s  s 


and  c  is  the  velocity  of  light. 

The  mean  time  delay  At^  of  the  portion  of  the  light  pulse  radiated  into 
the  n-th  element  of  solid  angle  is  given  by 


2  2  2 
(x  -  x  )  +  (y  -y  )  +  (z  -  z  ) 

s  n  Js  Jn  s  n 


+  t  -  t 
n  o 


The  energy  arriving  at  the  satellite  from  the  element  of  cloud 
top  area  at  z^,  y^,  x^,  is  proportional  to  Tn  and  is  treated  here  as  arriving 
as  a  pulse  at  At^.  Normalization  is  done  in  terms  of  the  direct  energy 
the  satellite  would  receive  in  the  absence  of  clouds. 

When  all  N  points  have  been  treated  in  this  way  the  accumulated 
normalized  energies  in  the  time  boxes  are  plotted  sequentially  from  0  to 
299  nsec.  Such  plots  are  interpreted  as  representing  the  power  as  a  function 
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of  time  that  would  arrive  at  the  satellite  if  a  delta  function  light  source 

were  fired  at  the  location  and  the  time  of  the  cloud  photographs. 

Using  the  same  values  of  x  ,  y  ,  z  ,  At  ,  and  T  similar 
6  n  Jn  n  n  n 

power  vs  time  curves  can  be  calculated  for  satellites  in  other  directions 

than  that  of  the  sun.  (These  directions  involve  new  values  of  x  ,  y  ,  z  .) 

s  s  s 

Some  additional  assumptions  are  involved  when  this  is  done  but  they  are 
probably  not  as  great  as  those  already  made. 
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SECTION  3 


CERTAIN  THEORIES  AND  APPROXIMATIONS 


3.  1  TRANSPORT  MEAN -FREE  PATH 

Reasonable  analytic  statements  about  overcast  transmission 
and  associated  time  delays  can  be  made  using  the  concept  of  the  transport 
mean-free  path.  This  was  first  introduced  into  optical  computations  by 
Sigmund  Fritz  (Ref.  3)  who  adapted  the  idea  from  work  on  neutron  dif¬ 
fusion  done  by  Victor  Weisskopf  (Ref.  4).  This  powerful  computational 
tool  is  seldom  used  in  optical  calculations.  It  is  applied  here  to  the 
discussions  of  two-stream  theory  and  diffusion  theory  outlined  below. 
Edward  Bucher  uses  the  concept  of  transport  mean-free  path  and  discusses 
some  of  its  limitations  in  his  paper  on  Monte  Carlo  simulations  of  pulse 
propagation  through  clouds  (Ref.  5). 

A  simple  overcast  may  be  represented  as  a  layer  of  great  lateral 
extent  and  of  thickness  B  wherein  there  are  uniformally  distributed  scatter¬ 
ing  centers  resulting  in  an  average  scattering  coefficient,  a,  and  phase 
function  for  scattering  <&  (  4>  ).  Corresponding  to  this  realistic  represen¬ 
tation  there  is  a  pseudo-overcast  having  the  same  thickness  but  made  up 
of  isotropic  scattering  centers  of  size  and  number  density  such  that  in  this 
pseudo-medium  the  scattering  coefficient  £  is  given  by 


where 


cos  <t>  = 


2  =  ct(1  -cos  <t> ) 

£*(  +  )  cos  <f>  sin  4>  d 


/’ 


<$>  (  <(>  )  sin  <f>  d  <£ 


3.  1 


3.  2 
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that  is  cos  <t>  is  the  average  cosine  of  the  scattering  angle  relative  to  the 
direction  of  the  incident  radiation.  In  the  "real"  medium  the  mean -free 
path,  l  is  given  by 


1=  1/a  3.3 

In  the  pseudo-medium,  the  value  L  is  called  the  transport  mean-free 
path  and  is  given  by 

L  =  -5-  =  l  ( 1  -  cos  <t>  )  1  3.  4 

In  some  diffusion  theory  calculations  the  mean -free  path  and  its  associated 
phase  function  may  be  replaced  by  the  transport  mean -free  path  and  its 
associated,  single  valued  phase  function  which  is  equal  to  1/4  n.  This  can 
result  in  a  major  reduction  in  the  complexity  of  calculations.  Bucher  (Ref. 
5)  shows  that  when  the  ratio  of  overcast  thickness  B  to  transport  mean- 
free  path  L,  B/L,  is  larger  than  2  or  3  then  the  isotropic  scattering  approx¬ 
imation  may  be  used  in  diffusion  calculations  without  serious  compromise 
of  accuracy. 

In  the  work  reported  here  measurements  are  made  and  inter- 
2 

preted  in  terms  of  B/ L  or  B  / L.  For  instance,  measured  overcast 

transmission  is  used  to  calculate  B/L.  Measured  values  of  oxygen 

2 

absorption  within  overcasts  are  used  to  calculate  values  of  B  /  L.  The 
work  done  here  does  not  deal  with  £  ,  the  true  mean-free  path  in  the  over¬ 
cast  or  with  (ttST  <f>  ,  the  average  cosine  of  the  scattering  angle. 

3.  2  TWO-STREAM  THEORY 

Two -stream  theory  may  be  used  to  develop  expressions  for  the 
transmission  and  the  reflection  of  a  uniform  overcast  of  great  extent.  When 
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it  is  assumed  that  there  is  no  absorption  within  the  overcast  and  no  re¬ 
flection  from  the  ground  below  the  overcast,  equations  such  as  those  given 
by  Schuster  (Ref.  6)  and  Dietzius  (Ref.  7)  reduce  to  those  given  by  Mecke 
(Ref.  8).  Mecke 's  methods  of  calculations  have  been  applied  to  the  case  of 
isotropic  scattering  and,  for  this  simplification,  equations  have  been  de¬ 
veloped  for  the  radiance  as  a  function  of  angle  of  view  of  the  bottom  and  of 
the  top  of  an  overcast  for  any  angle  of  incidence  of  sunlight  (Ref.  9  &  10). 
Hiese  results  are  summarized  here. 

For  a  uniform  overcast  of  great  extent  made  up  of  isotropic 
scattering  centers  with  no  absorption  within  the  overcast  and  no  ground 
reflection,  the  transmission  of  the  overcast,  T,  when  illuminated  by  diffuse 
radiant  power  is 


T  =  1  +  b/l 


3.  5 


where  B  is  the  overcast  thickness  and  L  is  its  transport  mean -free  path. 
For  diffuse  illumination  the  reflection,  R,  is 


R 


B/L 

1  +  b/l 


3.  6 


When  parallel  radiant  power  is  incident  on  the  overcast  at  an 
angle  z  from  the  normal  to  the  overcast  the  diffuse  reflection  is  given  by 


R  ■ 
z 


2—  -  (  2  cos  z  -  1  )  |  1  -  exp  (  --|j  sec  z)  ^ 


2  (  1  +  B/L) 


3.7 


23 


and  the  diffuse  transmission  is 


2  +  (  2  cos  z  -  1)  |  1  -  exp  (  sec  z)  | 


2(1  +T> 


When  parallel  radiant  power, H  ,  is  incident  on  the  overcast 

o 

at  angle  z  from  the  normal,  the  radiance  of  the  top  of  the  overcast,  N 

tz 

is  independent  of  angle  of  view  and  equal  to 

^  cos  z  |  2-^-  -  (  2  cos  z  -  1)  |  1  -  exp(-^r-  sec  z)jj  3.9 

tZ  tt  \  2(  1  +  B/L)  / 

The  radiance  of  the  base  of  the  overcast,  ,  is  dependent  on  the  angle 

of  view,  Q,  measured  from  the  downward  normal  to  the  base.  It  is  given 

by 

3H  cos  z  (2  cos  6  +  1)  i  2  +  (2  cos  z  -l)(l-exp  [  ~  sec  z  ]  )  1 

bz0  =  — - 1 - - h - L  o  , 

14  t t  (  1  +  B/L) 


When(B/L)sec  z  »  oo,  that  is,  for  overcasts  through  which 
one  cannot  see  the  sun 
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3.  3  THE  REDUCED  ATMOSPHERE 

Considerable  simplification  in  the  working  equations  and  cal¬ 
culations  used  in  describing  scattering  and  absorption  in  clouds  result 
from  transferring  the  problem  from  the  real  to  the  reduced  atmosphere. 
The  reduced  atmosphere  is  one  of  constant  density  equal  to  that  at  sea  level 
and  extending  from  the  sea  level  to  an  altitude  H  called  the  scale  height. 
Above  this  height  the  atmospheric  density  is  zero.  The  same  atmospheric 
mass  above  unit  area  at  sea  level  is  involved  in  the  real  as  is  the  reduced 
atmosphere.  Within  the  reduced  atmosphere  the  molecular  scattering 
coefficient  and  the  oxygen  absorption  coefficient  are  constant. 

It  is  assumed  here  that  in  the  true  atmosphere  the  air  density, 
p  (h),  is  an  exponential  function  of  height,  h,  obeying  the  law 


P(h) 


p  (0)  e 


-h/H 


3.  15 


where  P  (0)  is  the  atmospheric  density  at  sea  level  (h  =  0 )  and  H  is  the 

scale  height.  Unless  stated  otherwise  H  is  taken  here  as  equal  to  8  km 

_  3  _  3 

and  P  (0)  is  taken  to  be  1.  29  x  10  g  cm  .  With  this  model  for  the  real 
atmosphere  the  reduced  atmosphere,  having  vertical  coordinate  Z,  has 
constant  density  equal  to  P  (0)  from  Z  =  0  to  Z  =  H  and  zero  density  for 
Z  >  H.  That  is 

P  (Z)  =  P( 0)  for  0  <  Z  <  H 

and  p  (Z)  =  0  for  Z  >  H  3.  16 
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Values  of  Z  and  h  are  related  as  follows 


and 


Z 


H  (1  -  e 


-h/H 


dZ  -h/H 
~6h  =  6 

h  =  H  ln  (  TT^zT* 
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If  a  cloud  in  the  real  atmosphere  has  its  base  at  altitude  h  and 
its  top  at  altitude  h  +  B  then  in  the  reduced  atmsophere  the  thickness  of 
the  cloud  j3  is  given  by 


/ 3  = 


H  e 


-h/H 


(  1 


-B/H. 
e  ) 
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If  a  cloud  in  the  reduced  atmosphere  has  its  base  at  altitude  Z  and  its  top 
at  altitude  Z  +  0  then  the  thickness  in  the  real  atmosphere,  B,  is  given  by 


B 


H  ln 


(  H  -  Z  ) 
(11-7-0) 
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3.4  DIFFUSION  THEORY  AND  MONTE  CARLO  CALCULATIONS 

The  probability  as  a  function  of  time,  p(t),  that  a  photon  in¬ 
cident  on  the  top  of  an  overcast  of  great  extent  at  time  t  =  0  will  leave 
its  bottom  at  time  t  may  be  stated  analytically  by  use  of  stochastic  methods 
(Ref.  11)  or  diffusion  theory  (Ref.  12),  or  may  be  calculated  by  Monte 
Carlo  simulation  (Ref.  5).  If  the  scattering  processes  in  the  overcast  arc 
modeled  as  isotropic  and  of  uniform  distribution  a  single  value  of  transport 
mean-free  path  may  be  assigned.  When  this  is  done,  the  analytic  statements 


reduce  to  forms  which  can  be  solved  easily  and  the  Monte  Carlo  simulation 
can  be  carried  out  with  greatly  reduced  effort.  The  diffusion  theory 
equations  derived  by  J.S.  DeGroot,  R.M.  Frank  and  H.S.  Stewart  (Ref.  12) 
are  given  here  with  a  few  illustrative  calculations  of  p(t^) .  Here  the  time 
base  is  changed  so  that  t^  =  0  at  the  earliest  time  that  a  photon  could  exit 
from  the  bottom  of  the  cloud  .  Then  t^  =  t-B/c  j  1  -  m  (  1  -  cos  z)  |  . 

Monte  Carlo  results  obtained  by  E.A.  Bucher  (Ref.  4)  are  compared  to  the 
analytic  results  and  certain  useful  relationships  developed  by  Bucher  are 
simplified  by  use  of  two-stream  theory  and  presented  here. 


The  equation  for  p  (t^)  At  developed  in  Reference  12  is  given 

here  with  some  changes  of  symbols. 

_  1 
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where 
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m 
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+  B  ,  » 

*  ~  |  l-m(l  -cos  z)  { 

=  m(cos  z  +  0,  71) 

=  (l-2n)-.71  m 

=  cloud  thickness,  km 
=  velocity  of  light 
=  L/B 

=  transport  mean-free  path,  km 
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Figure  3. 1  presents  plots  of  p  (t  )  vs  t  for  a  hypothetical  overcast  when 

a  _o 

the  value  of  B  used  was  1.  56  km  and  the  value  of  m  was  1. 92  x  10  .  The 

dotted  curve  is  calculated  using  diffusion  theory  and  the  solid  curve  is  from 
Figure  6  of  Bucher  (Ref.  5).  Figure  3.  2  is  similar  to  Figure  3.  1  except 
here  B  =  4. 14  km  and  m  =  0.0725.  Bucher  used  a  phase  function  for  scatter 
ing  calculated  by  Mie  theory  for  A  =  0.  7/u  and  for  a  distribution  of  water 
droplet  radii  between  3.  8y  and  4.  8/u.  The  minor  differences  between  the 
Monte  Carlo  curves  and  the  diffusion  theory  curves  appearing  in  Figures  3.  1 
and  3.  2  arise  from  the  difference  in  scattering  models.  In  Figure  3.  3 
the  scattering  phase  function  used  by  Bucher  (Mie)  is  compared  to  that  used 
in  the  diffusion  solution  (isotropic).  In  this  figure  the  more  conventional 
phase  function  for  scattering  <t>  (  4> )  has  been  multiplied  by  sin  <f>  .  The 
effect  of  the  strong  forward  scattering  in  the  Mie  theory  compared  to  the  iso 
tropic  scattering  of  the  diffusion  theory  model  seems  very  minor  for  the 
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Figure  3.  3.  Comparison  of  the  Mie  phase  function  to 

the  isotropic  phase  function  used  in  diffusion 
theory. 

relatively  thick  clouds  of  Figures  3. 1  and  3.  2  (diffuse  transmissions 
0.  25  and  0.  11).  However,  the  effect  is  quite  pronounced  in  an  optically 
thin  cloud  such  as  that  of  Figure  3.  4  where  the  diffuse  transmission  is 
42  percent.  Although  the  Monte  Carlo  results  are  more  nearly  representa¬ 
tive  of  real  overcasts,  diffusion  calculations  are  used  here  because  of  their 
low  cost  in  time  and  money. 
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Figure  3.4.  Time  dependent  probability  of  delayed  transit  for  a  hypo¬ 
thetical  optically  thin  overcast. 
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In  equation  3.  20,  which  describes  the  time  dependent  probability 

of  delayed  transit  of  a  photon  incident  on  an  overcast,  there  is  a  term, 

e  6  ,  which  accounts  for  absorption  processes  which  may  occur  within  the 

cloud.  The  curves  of  Figures  3.  1,  3.  2  and  3.  4  were  calculated  with 

g(t)  =  0,  that  is,  with  no  absorption.  In  the  case  of  the  oxygen  absorption 

treated  here  g(t)  =  b  {  /3/B  (c  t)  +  Y  }  since  t  =  nAt.  For  the  data  re- 

-1/2 

ported  here  b  =  0.  51  km  .  For  the  cloud  of  Figure  3.  1  B  =  1. 56  km 
(/3  =  1.  33  km  in  the  reduced  atmosphere)  and  if  the  zenith  angle  of  the  sun  is 
O5  and  that  of  the  tine  of  sight  603  then  the  total  path  through  oxygen  above 
and  below  the  cloud,  Y  is  7.  15  km.  The  term  g(t)  includes  Y,  the  external 
path  through  oxygen  because  the  fraction  of  the  total  absorption  which  occurs 
within  the  cloud  is  dependent  on  the  fraction  absorbed  external  to  the  cloud. 
Figure  3.  5  presents  a  plot  of  p(t^)  versus  t  for  a  hypothetical  cloud  including 
absorption.  Note  that  absorption  reduces  the  mean  time  of  delay  from 
14.4psec  (without  absorption)  to  12.  7  p sec. 

Equation  12  from  Bucher's  paper  (Ref.  5)  is  rewritten  here  as 

.  0.94 

cT  =  0.  62  B  (— ) 
d  m 

or  t.  =  —  -  I  c  =  0.3km/psec  3.21 

d  0. 94 

m 


Equation  17  of  the  same  reference  is  rewritten 
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A 

o 

C  +  1/m 


3.  22 
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Figure  3.  5.  Comparison  of  the  probability  of  delayed  transit  with  and 
without  C>2  absorption. 
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Bucher  notes  quite  a  pedigree  for  this  expression  which  he  atributes  to 

Van  de  Hulst  (Ref.  13)  and  Danielson,  Moore  and  Van  de  Hulst  (Ref.  14). 

The  value  of  C  is  given  as  1.42.  Values  of  A  as  a  function  of  z  ,  A  (  z), 

o  o 

are  given  in  Table  3. 1  from  Danielson  et  al  (Ref.  14).  Included  in  Table  3.  1 
are  values  of  Ao{z)  calculated  from  two-stream  results,  see  equation  3.  12, 
and  normalized  to  Danielson  et  al  at  cos  z  =  1.  Using  this  approximation 
and  equations  3.  21  and  3.  22  one  may  write  the  following  sequence  of  equations. 

-  1-13  (cos  z  +  1/2) 

'  1.42  +  — 

m 


_L_ 

m 


1.13  (cos  z  +  1/2) 


-  1.42 


Td  =  t  im.  1,42) 


0.  94 


Jd_=  2.36£OS-l-+..y*  -  1.26)°- 94 


B 


or,  in  terms  of  time  measured  so  that  incidence  on  top  of  cloud  occurs 
at  r  =  0 


2.36  (.™L?.U/2  -  !.  26,°-  94  +-L 
b  r  c 


3.  23 


This  hybrid  equation  plays  an  important  part  in  data  reduction. 


Table  3. 1.  Comparison  of  values  of  Aq(z)  reported  by  Danielson  et  al 
and  calculated  by  two -stream  theory. 

A  (z) 
o 

From  Reference  14 

1.  69 
1.  59 
1.  38 
1.16 
0.  92 
0.  62 

! 

0.52  -  0.62 


3.  5  LOW  RESOLUTION  02  -  BAND  ABSORPTION 

o 

In  the  work  reported  here  the  0„  absorption  in  a  10  A  band 
centered  on  .  7  606p  is  measured  for  the  path  sun  through  overcast  to  point 
observed  on  its  base  and  on  to  detector.  The  results  so  obtained  are  ex¬ 
pressed  in  terms  of  multiples  of  the  scale  height  of  the  atmosphere.  This 
choice  of  units  comes  from  the  method  of  calibration  wherin  the  absorption 
to  the  sun  is  measured  on  a  clear  day  from  sunrise  to  noon  or  from  noon 
to  sunset.  This  is  a  variation  on  the  Langley  method.  The  results  of  the 


1  fid 

—  (cos  z  +  1/2) 

1 .  D 

1. 69 
1. 58 
1.  35 
1.  13 
0.90 
0.  68 
0.  62 
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1 

0.9 
0.  7 
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0.3 
0.1 
0.05 


calibration  fit  an  expression  for  oxygen  transmission. 


T  ,  of  the  form 

n 


-b  yjH  sec  z 
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where  H  =  scale  height  of  atmosphere  =  8  km 
z  =  zenith  angle  of  sun 


The  square-root  relationship  between  In  Tn  and  the  sea  level  equivalent 

z 

path  length  through  O  ,  H  sec  z,  is  not  surprising.  It  is  discussed  briefly 

z 

in  an  excellent  outline  of  molecular  band  absorption  theory  which  is  given 
in  McClatchey  et  al  (Ref.  15)  along  with  references  which  treat  the  subject 
in  depth. 

In  the  reduced  atmosphere  the  oxygen  density  is  constant  and 
path  length  through  atmospheric  oxygen  is  equal  to  path  length  per  se. 
Elapsed  time  and  path  length  through  atmospheric  oxygen  are  related  through 
c  the  velocity  of  light.  Equation  3.  24  may  be  written 


-b  }j  26.  67  c  sec  z 


There  is  a  mean  transit  time  for  the  diffusion  of  light  through  an  over¬ 
cast  given  by 


3.  25 
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r 
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where  p  (t)  is  the  probability  that  a  photon  will  penetrate  the  overcast  in 
time  t.  In  the  reduced  atmosphere  the  mean  time  of  transit,  r  ,  is 

“  =  t  /3/B 


where  B  is  the  true  cloud  thickness  and  3  is  that  thickness  in  the  reduced 
atmosphere.  The  probability  function  in  the  reduced  atmosphere,  p  (  r  ), 
is  related  to  that  of  the  real  atmosphere  as  follows 


p  (r)  *  p(t)  /3/B 

/q00  P  <t>t  dT 

f<T p  {t)  dT 


The  transmission  of  sunlight  through  the  overcast,  T,  is 


T  =  p(t)  dr 


3.27 


The  transmission  of  oxygen  in  the  total  path  sun  through  overcast  to  detector, 

-b  ifi 


T  ,  is  written  as 
°2 


Y  +  cr 


Or 


jfc0p<T>d’ 
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where,  in  the  reduced  atmosphere,  Y  is  the  path  length  of  sunlight  from  the 
top  of  the  atmosphere  to  the  top  of  the  cloud  plus  the  path  length  from  the 
observed  point  on  the  bottom  of  the  cloud  to  the  detector. 


1 
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By  the  definition  of  *  ,  equation  3.2  6,  the  term  in  parentheses  goes  to 
zero  and 


O, 


~  e 
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Y  +  cT 


3.  36 


Equation  3.  36  is  the  heart  of  the  work  reported  here  because  it  relates 
measured  oxygen  transmission  to  mean  transit  time  through  the  overcast. 
It  is  an  approximation  based  on  the  conditions  of  equations  3.  30  and  3.  34. 

It  is  a  good  approximation  if  for  all  times  for  which  p  (r)  is  significant  the 
following  conditions  hold 


and 
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The  photographs  taken  in  the  work  reported  here  were  re¬ 
duced  with  equivalent  resolution  of  0.43  degrees  in  space.  It  is  probable 
that  most  of  the  data  taken  meets  the  above  conditions.  This  is  due  in  part 
to  the  small  solid  angle  represented  by  each  pixel.  The  conditions  may  not 
be  met  at  the  edges  of  clouds  or  in  the  presence  of  complex  overcasts.  The 
risk  of  violating  the  assumptions  increases  with  field  of  view  of  the  resolution 
element. 


SECTION  4 


RELATIONSHIP  USED  IN  DATA  REDUCTION 

4.  1  HEIGHT  OF  CLOUD  BASE 

The  heights  of  the  bases  of  the  clouds  were  determined  by 
analysis  of  stereo  pictures.  The  stereo  cameras  were  located  at  the 
primary  data  station,  A,  and  the  secondary  station,  B.  Station  A  was 
at  location  UCS2  at  Kennedy  Space  Center.  The  coordinates  of  A  were: 

Latitude  28°  31'  31.27" 

Longitude  80°  37'  20.61" 

The  separation  of  B  from  A  was  3215.5  feet  or  980  meters.  The  line 
from  A  to  B  was  rotated  clockwise  34°  01'  from  true  north.  The  stereo 
cameras  were  positioned  at  station  A  and  B  such  that  a  line  from  the  top 
to  the  bottom  of  the  image  format  was  parallel  to  the  line  from  station  A 
to  B  with  both  formats  oriented  the  same. 

Two  stereo  camera  pairs  were  used.  The  first  involved  Flight 
Research  cameras  with  18.  5mm  Angenieux  Lenses.  These  operated 
satisfactorily  but  their  limited  field-of-view  made  them  almost  useless  on 
the  many  occasions  when  cloud  bases  were  low.  The  second  stereo  pair 
used  Automax  cameras  with  Fairchild  160°  f/1.5  lenses  for  35mm  format. 
These  gave  satisfactory  coverage  for  almost  all  cloud  base  heights  en¬ 
countered.  Relationships  used  in  reduction  of  data  from  these  cameras 
are  discussed  in  the  following  paragraphs. 

Figure  4.1  presents  the  manufacturers'  relationship  between 
the  off-axis  distance  in  inches,  d,  in  the  image  plane  (called  the  image 
height  by  the  manufacturer)  versus  the  external  off-axis  angle,  z,  in 
degrees  for  the  Fairchild  f/1.5  lens.  An  empirical  relationship  matching 
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Image  Height,  d  (Inches) 


the  curve  of  Figure  4.1  is: 


9  4  2  “id 

z  =  1.373  x  10  d  +  0.12(10  ‘  -1)  4.1 

In  use  the  cameras  were  mounted  with  optical  axis  pointed  to  the  zenith 
and,  under  these  conditions,  the  value  of  z  assigned  to  a  particular  point 
was  the  zenith  angle  of  that  point.  A  photograph  of  an  Automax  camera 
with  Fairchild  lens  attached  is  shown  in  Figure  4.  2.  A  picture  of  the  camera 
system  mounted  in  a  weatherproof  enclosure  is  given  in  Appendix  A.  The 
camera  and  enclosure  are  supported  by  a  pedestal  which  provided  adjustments 
at  its  base  to  orient  the  optical  axis  of  the  camera  system  vertical.  Such 
alignment  was  achieved  using  a  special  lens  cap  to  which  a  sensitive  level 
dome  was  attached. 

During  data  runs  an  operator  at  station  B  was  in  telephonic  con¬ 
tact  with  station  A  and  took  synchronous  pictures  with  those  taken  at  station 
A.  Proper  film  labeling  allowed  synchronized  stereo  pictures  to  be  matched 
after  processing.  Figures  4.  3  and  4.4  are  black  and  white  positives  made 
from  a  matched  pair  of  35mm  color  positive  transparencies  from  station 
A  and  station  B  respectively.  Particular  features  were  identified  on  both 
prints  and  numbered.  The  original  35mm  transparencies  were  viewed  with 
a  Richardson  film  reader  and  the  magnified  film  coordinates  of  the  numbered 
features  were  read  and  recorded.  The  images  were  rotated  in  the  reader 
using  the  film  edge  as  control  to  bring  each  image  of  the  stereo  pair  in  line 
with  the  camera  axis  defined  by  the  line  from  station  A  to  station  B.  Dis¬ 
placements  on  the  film  were  read  with  the  center  of  the  sky  image  as  the 
origin  of  an  x,y  coordinate  system.  With  knowledge  of  the  magnification  of 
the  Richardson,  image  height  values,  h,  could  be  calculated  and  used  in 


* 


Equation  4.  1  to  obtain  the  zenith  angle,  z,  or  a  particular  feature  as 
observed  by  the  camera.  The  azimuth  angle,  a  ,  of  the  feature  was  easily 
calculated  referenced  to  the  station  A-station  B  axis  from  the  feature 


coordinates.  Figure  4.  5  shows  the  relationships  between  a  ,  a  ,  z 

A  BA 

and  z  for  a  particular  feature  location,  F.  If  R  is  the  separation  of  A 
B 

and  B  then  h  the  height  of  feature,  F,  is  given  by 

R  cos  z„ 


h  = 


cos  a  -  sin  a „  cot  a  . 
B  BA 


4.  2 


and 
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tan  z  cos  a-  -  tan  zA  cos  o, 
B  B  A  A 


4.  3 


where  equation  4.  2  is  the  more  accurate  of  the  two  in  the  two  quadrants 
defined  by  a  between  45°  and  135°  and  a  between  225°  and  315°.  Equation 
4.  3  is  better  in  the  other  two  quadrants.  Values  of  the  readings  from  the 
Richardson  appropriate  to  Figures  4.  3  and  4.  4  are  presented  in  Table  4.  1 
with  the  calculated  values  of  h  using  both  equation  4.  2  and  4.  3. 

4.  2  CALIBRATION  OF  ALL  SKY  RADIOMETER  AND  OF  ALL  SKY  CAMERAS 
4.  2.  1  Geometric  Relationships 

The  basic  tools  of  this  experiment  record  data  reflected  from  a 
spheric  mirror  surface.  The  all  sky  radiometer  and  the  two  camera  systems, 
the  on-band/ off-band  and  the  red-blue  cameras,  are  included  in  this  state¬ 
ment.  Figure  4.  6  shows  such  systems.  The  camera  shown  in  Figure  4.  6a 
photographs  the  virtual  image  of  the  hemisphere  of  the  sky  which  is  formed 
somewhat  below  the  surface  of  the  reflecting  sphere.  In  Figure  4.  6b  a 
similar  virtual  image  is  re -imaged  by  an  objective  lens  onto  a  relay  lens 
assembly  which,  in  turn,  images  the  objective  lens  onto  a  radiation  detector. 
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Table  4.  1.  Cloud  base  heights  for  8/1/77  at  1715  hours  calculated  for 

the  numbered  features  in  the  stereo  pair  of  Figures  4.  4  and  4.  5 


Feature 

No. 


Richardson  Readings 


Feature  Height 


Station  A 

Station  B 

by  Eqn  1 

X 

y 

X 

y 

meters 

meters 

706 

-431 

729 

-254 

4800 

4715 

450 

19 

434 

241 

4705 

4714 

241 

-273 

234 

-5 

4285 

4056 

270 

174 

262 

401 

4539 

4504 

-53 

-366 

-60 

-126 

3908 

4441 

-254 

+347 

-238 

+577 

3955 

3925 

Arbitrary  Units 


COORDINATE  SYSTEM 


A  ft  8 

Figure  4.  5.  Geometric  relationship  of  feature  f  as  seen  from  Stations 
A  and  B  for  cloud  height  determination. 


P 
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Detector 
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Reflecting 

Sphere 
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All  Sky  Radiometer 


Figure  4.  6.  All  sky  detectors  using  reflections  from  a  spheric  surface 


With  reference  to  the  camera -sphere  system  of  Figure  4.  6a 
the  relationship  between  the  external  angle  of  the  camera  0  and  the  external 
angle  of  the  camera  sphere  system,  z,  is  shown  in  Figure  4,  7.  In  that 
figure  the  optical  axis  of  the  camera  is  directed  vertically  downward 
toward  the  center  of  the  sphere  at  C.  The  radius  of  the  sphere  is  R  and 
the  distance  from  the  sphere  surface  to  camera  lens  is  D.  A  line  from  the 
lens  center  with  external  angle  0  intercepts  the  spheric  surface  at  point  P 
and  is  reflected  at  angle  z  from  the  vertical.  The  line  from  P  to  the  center 
of  the  sphere  is  at  angle  a  to  the  line  from  sphere  center  to  center  of  the 
camera  lens  (the  vertical).  The  horizontal  distance,  y,  of  the  point  P  from 
the  line  between  lens  center  and  sphere  center  may  be  written  as 

y  =  R  sin  a 


or 


y  =  [D  +  R  (1  -  cos  a)]  tan  0 


and  so 


0  =  tan 


1 


R  sin  a _ 

D  +  R  (1  -  cos  a) 


By  inspection  of  Figure  4.  7  one  may  write 


4.  4 


z 


2  o+0  =2  a  +  tan 


R  sin  a 

D  +  R  ( 1  -  cos  a) 


4.5 


As  the  ratio  of  R/D  goes  to  0  the  angle  0  goes  to  0  and  a  =  z/2.  That  is, 
when  z  =  90°  and  the  line  from  the  lens  to  P  is  reflected  to  the  horizon, 

a  •  45°.  in  practice  R/D  is  finite.  If  it  is  equal  to  0.  1  then  the  value  of 

h 

a  such  that  the  line  from  the  lens  is  reflected  to  the  horizon  is  a,  =  43. 1° 
h  h 

and  if  R/D  =  0.  2,  a  =  41.4°.  Figure  4.  8  shows  the  relationship  between 
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Ratio  r/rgo 


the  external  angle,  z,  to  the  image  height  for  the  ideal  camera -sphere 

system  where  R/D  - *0.  What  is  plotted  in  Figure  4.8  is  the  ratio 

of  the  image  height  at  external  angle  z  to  the  image  height  at  z  =  90°. 
For  comparison  Table  4.  2  presents  the  deviation  of  this  ratio  from  the 
ideal  case  when  R/D  =  0.1  and  R/D  =  0.2.  It  is  apparent  that  minor 
errors  are  introduced  by  treating  the  system  with  R/D  <0.2  as 
ideal  (i.  e.  R/D  =  0). 

Table  4.2.  Magnitude  of  errors  introduced  by  treating  finite  camera- 
mirror  separation  as  infinite. 


External 

Angle 

Degrees 

Fractional  Deviation 

R/D  =  0.1 

R/D  =  0.2 

5 

1.011 

1.020 

15 

1.012 

1.020 

30 

1.011 

1.017 

55 

1.008 

1.013 

75 

1.004 

1.006 

In  the  following  we  discuss  the  relationship  between  area  in 
the  focal  plane  of  a  camera-sphere  system  and  external  solid  angle.  The 
approximation  z  =  2 a  ,  which  was  discussed  above,  is  used.  With  ref¬ 
erence  to  Figure  4.7,  a  point  in  the  focal  plane  of  the  camera  can  be 
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located  Ln  the  polar  coordinate  system  p  0.  In  this  case 

P  =  M  R  8in  a  4.6 

where  M  is  the  magnification  of  the  camera  as  used  and  M  a  F/D  The 
angle  0  is  about  the  optical  axis  and  is  measured  from  any  convenient 
fiducial  (usually  north).  As  shown  in  Figure  4.  9  an  element  of  area  in 
the  focal  plane,  da,  is 

da  =  p  dp  d0 
2  2 

da  =  M  R  sin  a  cos  a  da  do  4.7 

Referring  again  to  Figure  4.  9,  the  solid  angle 

dfi  =  —  sin  z  dz  d0  4.8 

but  z  =  2  a 

or  d  ft  =  s  in  2a  da  d  6  4.9 

Since  sin2»  =  2  sin#  cos  a 

4  10 

d ft  =  2  sin  a  cos  a  do  d0 

2 

or  dft  =  — « — 2~  da  ,  steradians  4.11 

M  R 

This  exercise  has  demonstrated  that  for  a  particular  camera -sphere  system 
there  is  a  single  valued  coefficient  relating  area  in  the  camera  focal  plane 
to  the  associated  external  solid  angle.  This  is  an  important  concept. 

For  example,  if  the  hemisphere  of  the  sky  is  photographed  by  such  a  camera- 
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Figure  4.9.  Relationship  between  an  element  of  area,  da,  in  the 

image  plane  (upper  drawing)  to  an  element  of  solid  angle, 
dO  ,  in  object  space  (lower  drawing). 


sphere  system  and  if  30  percent  of  the  exposed  circular  area  represents 
clear  sky  (as  opposed  to  cloudy)  then  there  is  30  percent  of  the  hemisphere 
of  the  sky  for  which  there  is  cloud  free  line -of-  sight. 

4.  2.  2  Camera  System  Radiometric  Relationship 

If  the  exposed  circular  area  in  the  focal  plane  of  the  camera 

recording  system  as  shown  in  Figure  4.  8  is  divided  into  P  elements  of 

solid  angle  A  a.  these  will  correspond  to  P  elements  of  solid  angle  An 

th 

in  the  hemisphere  of  the  sky.  The  exposure  E  of  the  p  element  A  ap 
is  given  by 


■n  N  At  T 
P _ 

4  (f/n)  2 


4.  12 


th 

where  Npis  the  radiance  of  the  p  solid  angle  element  of  the  hemisphere 

of  the  sky,  A  t  is  the  exposure  time,  T  is  the  transmission  of  the  camera- 

sphere  system  and  (f/n)  is  the  camera  f-number  (ratio  of  focal  length  to 

aperture  diameter).  The  product  of  the  exposure  and  the  area  of  an  element 

E  A  a  is  the  energy  incident  on  the  area  during  the  exposure.  The  total 
P  P 

energy  involved  in  an  exposure  is 


E 

p 


E  Aa  =  ~— 
p  p  4 


At 


t  mV 


(f/n  r 


E 

p 


N  A  Q 
P  P 


4.  13 


If,  as  in  Figure  4.  6b,  the  focal  plane  of  the  camera  is  replaced 
by  a  relay  lens  of  transmission  Tr  forming  an  image  of  the  objective  lens 
onto  a  detector,  then  the  energy  incident  on  the  detector  in  a  time  A  t  is 


T  D  ED  A 
r  p  P 

unit  time  or 


ap  .  The  power  incident  on  the  detector,  P  ,  is  the  energy  per 


w 

T 


2  2 

T  T  R  M 

- £ -  In  An 

p  P  p 


4.  14 


(f/n)- 
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If  the  only  source  of  power  were  the  sun  then 

H  T  ,  the  local  solar  irradiance.  Here  T  is 
©  ©  o 

mission  in  the  path  to  the  sun  and  is  the  solar  irradiance  outside  the 
earthb  atmosphere.  These  equations  are  the  basis  of  the  camera  systems 
calibration  through  the  all -sky  radiometer. 

4.  2.  3  All-Sky  Radiometer  Calibration 

Calibration  of  the  all -sky  radiometer  is  carried  out  in  the  following 
way.  The  filter  condition  is  selected,  this  implies  that  either  no  filter  is 
used  and  the  spectral  response  is  characterized  by  the  bare  silicon  detector, 
or  the  red  filter  is  used.  The  red  filter,  which  matches  the  red  filter  used 

O 

in  the  red-blue  all  sky  camera,  has  a  100  A  band  pass  centered  on  0.  672m. 

Calibration  with  that  filter  is  described  here.  The  transmission  to  the  sun, 

T  ,  at  the  wavelengths  0.44m,  0.  5m,  0.  64m,  0.  88m  and  1.  68m  is  measured 

using  a  Sun  Photometer  made  by  Fred  Volz  of  AFGL.  From  a  plot  of  these 

transmission  readings,  a  value  for  T^o  the  transmission  to  the  sun  at  0.  672m 

is  obtained.  A  reading,  s^  +  sg,  of  the  output  of  the  all  sky  radiometer  is 

made  as  it  is  exposed  to  the  full  sky,  s  ,  plus  direct  sun,  s  .  Then  the 

s  © 

direct  sun  is  shaded  so  that  only  the  sky  irradiates  the  radiometer  and  the 
reading  of  sky  signal,  s  ,  is  made.  The  reading  S  which  would  be  obtained 

S  0 

by  the  radiometer  if  it  were  outside  the  earth's  atmosphere  and  illuminated 
by  the  sun  above  is  given  by 


y'  N  AO  would  equal 
p  P  P 

the  atmospheric  trans- 


S 

o 


(s  +  s  )  -  s 
OS  s 


Kn  H  A  T. 

O  £>  I 


4.  15 


where  Kq  is  the  calibration  coefficient  of  the  all  sky  radiometer  and  filter, 
H  ^  sPectral  irradiance  of  the  sun  outside  the  earth's  atmosphere  at 

the  central  wavelength  of  the  filter,  AX  is  the  band  pass  of  the  filter  and  T^ 
is  the  mean  transmission  of  the  filter. 
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The  signal  from  the  sun  plus  the  sky  (s  +  s  )  is  related  to  the 

0  S 

radiance  of  the  sky  and  the  apparent  radiance  of  the  sun's  disk  as  follows: 


s 


s  +  s 
©  s 


K 


o 


K 

o 


T.  £  N  A  Q  AX 
f  p  P  P 


4.  16 


Relations  4. 15  and  4. 16  do  not  depend  on  the  details  of  the  all  sky  radio¬ 
meter  design,  however  the  value  of  K  does.  In  any  situation,  overcast, 
broken  clouds  or  clear  the  observed  signal,  s,  is  related  to  the  detailed 
radiance  of  the  hemisphere  of  the  sky  by  equation  4.  16.  In  connection  with 
these  concepts  the  term  "coupling  coefficient"  has  been  defined  as  follows 

T*  £  N  A  Q  AX 

Coupling  Coefficient  =  — t - £ B - E -  =  — 

s  s 

o  o 

The  same  reasoning  applies  to  the  case  of  the  all -sky  radio¬ 
meter  used  without  a  filter.  The  results  obtained  are  given  in  terms  of 
coupling  coefficients  appropriate  to  bare  silicon.  The  transmission  to  the 
sun  is  somewhat  more  complicated.  It  is  defined  as  follows: 


T 

© 


/H  ,  G  ..  T  .  dX 
J  ©X  siX  ©X 

7H  .  G  ...  dX 
©X  siX 


4.  17 


where 


H 

©X  is  the  spectral  solar  irradiance  outside  the  atmosphere. 
^8iX  *s  t*le  relative  spectral  response  function  of  silicon. 

T0x  *s  t*ie  sPectr&l  atmospheric  transmission  to  the  sun. 


The  values  ofH  >  and  G  . .  are  found  in  handbooks.  The  values  of  T  .  are 
©A  siA  ©A 

taken  from  the  plot  made  from  Sun  Photometer  readings. 
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4.2.4  Red  Channel  Camera  System  Calibration 

The  objective  of  the  calibration  of  the  all-sky  camera  (camera- 
sphere  system)  is  to  relate  through  a  calibration  factor,  K^,  the  measured 
relative  exposure  E  in  an  element  of  area  A  a  to  the  radiance  of  the  sky 

r  r 

Npin  corresponding  element  of  solid  angle,  Afi^.  The  relationship  is 
derived,  however,  from  the  total  image  area  and  total  sky  solid  angle  and 
their  respective  elements.  The  calibration  is  made  by  comparing  results 
from  the  red  all-sky  camera  records  with  readings  from  the  calibrated  all 
sky  radiometer  for  the  same  time.  Since  both  systems  are  linear  in  their 
responses  and  since  an  identical  filter  is  used  in  the  one  as  in  the  other  the 
band  pass  and  transmission  of  the  filter  can  be  omitted  from  the  calculations 
of  calibration.  In  addition  the  exposure  time  of  the  camera  was  unchanged 
between  calibration  and  use  so  it  too  may  be  omitted  from  the  calculations. 
The  total  energy  involved  in  an  exposure  of  the  all  sky  camera  is  given  by 
equation  4.13  and  is  rewritten  here  with  the  above  comments  in  mind  as 

K,  £  E  da  =  ?  N  do  4*  18 

1  P  p  p  P  p  p 

and  so 

K.  E  =  N  4.19 

1  p  p 

The  value  of  is  found  by  the  following  procedure.  A  daytime  condition 
with  overcast  covering  the  sun  is  selected  (the  direct  solar  image  "over 
exposes"  the  film).  A  picture  is  made  using  the  red  filtere  •  all  sky  camera 
and  at  the  same  time  a  reading  is  made  with  the  red  filtered  all  sky  radio¬ 
meter.  For  this  case  Np  A  Op  is  the  same  for  the  two  instruments.  The 
calibration  signal  sc  fromthe  all  sky  radiometer  is  given  by  equation  4.16 
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which  is  rewritten  here  as 


s 

c 


A  0 
P 


and  the  calibration  of  the  all  sky  radiometer  gives 


4.  20 


S 

o 


K 

o 


H 

o 


It  follows  that 


N  Afl 
P  P 


4.21 


4.  22 


The  film  exposed  in  the  red  filtered  all  sky  camera  is  processed  and  read 
in  a  controlled  way  so  that  the  relative  values  Ep  are  determined.  The 
sum  ^Epdap  is  calculated  and  is  designated  here  as  EA  .  That  is,  for 
the  calibration  exposure 


EA  =  £ 
P 


E  Aa 
P  P 


and  from  equations  4.  18  and  4.  22 


K 


1 


s 

c 


EA 


H 

o 


S 

o 


4.  23 


4.  24 
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4.  2.  5  Red  Channel  Camera  Calibration  Results 


This  paragraph  outlines  a  calibration  of  the  filtered  all-sky 
radiometer  and  the  subsequent  calibration  of  the  filtered  all-sky  camera. 
The  all-sky  radiometer  was  calibrated  at  1034  on  August  2,  1977.  Average 
readings  were  as  follows 

sun  plus  sky  (s  +  s  )  =  77  6.  5  meter  units 

o  s 

sky  s  =  129.0  meter  units 

s 

(s+s)-s=s  =  647.  5  meter  units 

os  s  © 


Transmission  to  the  sun  was  measured  at  1013  and  at  1036  using  the  sun 
photometer.  The  resultant  transmission  curves  are  given  in  Figure  4.  10. 
A  value  of  transmission  at  0.  672p  was  taken  to  be  0.  76.  The  value  of  S 

o 

the  solar  signal  in  the  absence  of  the  atmosphere  is 

s 


S  =  ~ 
o  T 


=  852  meter  units 


3-2 

Now  H  for  X  =  672/u  is  1.  57  x  10  wm  ju 
oA 


or 


K 


2  -1  -1 

=  0.  525  (meter  units)  m  w  m 


4.  25 


If  the  reading,  s  ,  for  an  overcast  were  220  then  Z)  NrtA  fin  =  220/0.  525  = 

-2  P  ^  _2 

419  w  m  iu.  The  average  radiance  of  the  overcast  is  67  w  m  /j  sr  and 

the  average  transmission  of  the  overcast  is  13  percent.  Note  thatIZ  A  np  =  2 

P 

for  the  all  sky  radiometer. 


Calibration  of  the  all  sky  camera  with  the  red  filter  is 
discussed  in  this  paragraph.  At  1702  on  August  2,  1977,  the  sky  was 

generally  overcast.  A  red  all  sky  picture  was  taken  and  the  sumZ^En  Aan  = 

2  P  p  p 

was  found  to  be  equal  to  1. 015  x  10  .  At  the  time  of  the  photograph  the  all 

sky  radiometer  read  s  =  195.  From  the  paragraph  above.s  =  852  and 

3-2-1  0 

Hq  =  1. 57  x  10  w  m  n  .  Equation  4.  24  gives 

s’  H  q  ,,  —2  — 1 

K  =  _ £_  =  9.  85  x  10 watts  g  Mr _ 

*  '  relative  exposure  unit 


EA  S 


This  value  of  Kj  is  used  in  data  reduction. 
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4.  2.  6  Calibration  of  the  ON/ OFF  Band  Camera  System 

The  following  paragraphs  are  concerned  with  the  calibration 
of  the  filters  used  with  an  all-sky  camera  to  measure  the  oxygen  absorp¬ 
tion  associated  with  each  element  of  solid  angle  do  .  This  camera  system 
need  not  be  calibrated  in  absolute  terms  since  only  ratios  of  exposures  are 
used  in  data  reduction. 

For  each  experiment  two  pictures  were  taken  on  High  Speed 
Infrared  Film  with  the  all-sky  camera  assigned  to  oxygen  absorption  measure- 

o 

ments.  One  picture  was  taken  through  a  filter  having  13.  8  A  bandwidth 
at  half  maximum  centered  on  0.  7608p.  This  wavelength  interval  is  located 
within  an  atmospheric  oxygen  absorption  band  (the  A  band  extending  from 
its  band  head  at  0.  7593p  to  0.  778/j).  In  this  report  this  picture  is  referred 

o 

to  "on-band".  The  second  picture  was  taken  through  a  filter  having  33  A 
bandwidth  at  half  maximum  centered  on  0.7530p.  This  wavelength  interval 
is  free  of  any  absorption  by  gasses  normal  to  the  atmosphere.  This  picture 
is  referred  to  as  "off-band".  Figure  4.  11  presents  low  resolution  spectra 
of  the  transmission  of  the  atmosphere  for  three  lines  of  sight  having  zenith 
angles  0°,  60°,  and  70.  5°.  The  shape  of  the  envelope  of  any  curve  is  deter¬ 
mined  by  molecular  and  aerosol  scattering.  The  structure  of  the  curves  is. 
due  to  absorption  by  atmospheric  gasses.  The  location  of  the  peak  wave¬ 
length  of  both  the  on-band  filter  and  off-band  filter  are  shown.  Figure  4.  12 
is  a  higher  resolution  spectrum  of  the  A  band.  It  was  made  by  Curcio  and 
his  associates  at  Naval  Research  Laboratory  using  a  ten  mile  absorption 
path  across  Chesapeake  Bay.  The  location  of  the  center  of  the  on-band 
transmission  filter  is  shown.  In  this  figure  the  abscissa  represents  wave¬ 
length  and  the  ordinate  is  proportional  to  ( 1  -  T  )  where  T  is  the  trans- 

P  P 

mission  of  the  processed  spectrographic  plates. 


63 


OH  OH  OH  OlH  OH 


e  oxygen  absorption  band  of  interest 


It  was  anticipated  that  for  each  element  of  area  on  the  photo¬ 
graphic  films  the  ratio  of  the  exposure  through  the  on-band  filter  to  that 
through  the  off-band  filter  would  be  a  function  of  the  mean  quantity  of 
oxygen  in  the  path  camera  to  sun.  A  modification  of  the  Langley  method 
was  used  to  calibrate  this  ratio  using  known  paths  through  atmospheric 
oxygen.  The  filters  to  be  used  with  the  on-band/ off -band  cameras  were 
mounted  so  that  they  could  be  located  alternately  in  front  of  the  aperture  of 
a  Tektronix  J-16  Teleradiometer.  This,  in  turn,  viewed  a  white  plaque 
from  an  angle  45°  off  the  plaque  normal.  During  a  calibration  run  this 
assembly  was  positioned  so  that  the  sun  was  on  the  normal  to  the  plaque. 
Figure  4. 13  is  a  lay-out  sketch  of  the  assembly  and  Figure  4.14  is  a  photo¬ 
graph  of  the  instrument. 

The  equivalent  sea  level  path  length  through  atmospheric  oxygen  . 
L,  was  taken  to  be 


L  =  H  sec  z 


4.  26 


where  H 
z 


Scale  height  of  atmosphere 
Solar  zenith  angle 


The  scale  height  was  taken  to  be  eMual  to  8  km.  Calibration  ratios  taken  on 
3  August  1977  are  plotted  against  the  square  root  of  L  in  Figure  4. 15.  The 
equation  accounting  for  the  straight  line  through  the  points  in  this  figure  is 


or 


on  _  c  -0.51  Vl” 
off  6  e 


4.  27 


4.28 
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Figure  4.  13.  Schematic  diagram  of  calibration  equipment  for 
on-band/ off-band  camera  filters. 


Iters. 


Figure  4.  15.  Data  obtained  with  calibration  equipment  showing  the  logarithm  of  the 
on-band/ off-band  ratio  plotted  against  the  square  root  of  the  oxygen 
path  length. 


The  ratio  on/ off  is  used  to  indicate  the  ratio  of  the  teleradiometer  reading 
through  the  on-band  filter  to  that  through  the  off-band  filter.  Equations  4.27 
and  4.  28  are  used  in  data  reduction  to  account  for  the  ratios  of  calculated 
exposures  per  area  element  from  the  cn-band  and  off -band  filtered  all- 
sky  camera  photographs. 
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SECTION  5 


WORKING  EQUATIONS  FOR  DATA  ANALYSIS 

5.1  GENERAL  COMMENTS 

The  experiments  described  here  produce  a  limited  amount 
of  information  for  any  area  element  on  the  data  film  and  therefore  its 
corresponding  solid  angle  element  in  space.  For  each  element  the 
information  is  as  follows: 

a.  Whether  the  element  views  clear  space  or  a  cloud. 

b.  If  it  is  a  cloud,  a  value  for  the  height  of  its  base. 

c.  The  radiance  in  a  red  wavelength  band  of  the  portion 
of  the  cloud  viewed  by  the  element. 

d.  The  transmission  in  a  near  infrared  wavelength  band  of 
the  atmospheric  oxygen  in  the  mean  path  from  the  sun 
to  the  all-sky  camera  via  the  cloud  element. 

e.  The  azimuth  and  zenith  angle  to  the  sun  and  of  the  line 
of  sight  defined  by  the  element. 

There  is  such  an  assembly  of  information  for  each  of  the  127,000  elements 
of  area  on  the  film.  Data  reduction  relating  to  the  flow  of  power  in  the 
individual  elements  of  solid  angle  does  not  involve  interaction  of  information 
acquired  at  one  element  of  solid  angle  with  information  acquired  at  another. 
This  fact  reduces  the  complexity  of  the  bulk  of  the  data  reduction  procedure 
to  127,000  repetitions  of  rather  simple  calculations. 

It  seems  clear  that  the  radiance  of  an  area  on  a  cloud  is  closely 
related  to  the  efficiency  of  the  path  through  that  area  in  terms  of  getting 
light  from  the  sun  to  the  camera  station.  It  is  also  sensible  that  photon 
flight  times  involve  the  geometric  distances  from  the  camera  station  to 
the  area  viewed  and  from  the  corresponding  area  directly  above  on  the 
top  of  the  cloud  to  the  sun  in  addition  to  the  flight  time  will  in  the  cloud. 
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Finally,  absorption  of  light  by  oxygen  is  related  in  a  somewhat  compli¬ 
cated  way  to  path  lengths  in  the  atmosphere  within  and  free  of  clouds. 
These  complications  arise  from  the  altitude  dependence  of  oxygen 
pressure  and  the  square  root  law  of  absorption  by  oxygen.  There  may 
be  many  sets  of  simplifying  assumptions  that  allow  reasonable  calcula¬ 
tions  to  be  made.  The  set  used  here  are: 

a.  for  every  cloud  represented  by  a  data  element,  there 
is  a  corresponding  pseudo  cloud.  The  pseudo  cloud 
"exists"  in  the  reduced  atmosphere,  one  in  which  the 

pressure  of  O ^  is  constant  throughout.  The  observables 
of  the  pseudo  cloud  are  identical  to  the  real  cloud.  These 
observables  include  the  on-band/ off-band  ratio  and  the 
radiance  at  the  bottom  of  the  cloud.  The  characteristic 
of  the  pseudo  cloud  is  that  all  scattering  within  it  is 
isotropic. 

b.  the  pseudo  cloud  has  sufficient  horizontal  extent  that  edge 
effects  can  be  neglected;  that  is,  the  cloud  can  be  treated 
as  though  it  were  part  of  a  uniform  overcast  whose  proper¬ 
ties  are  defined  by  the  properties  of  the  element. 

c.  the  density  of  scattering  agents  within  the  pseudo  cloud 
in  the  reduced  atmosphere  is  constant  between  the  upper 
and  the  lower  boundaries  of  the  cloud. 

d.  the  size  distribution  of  the  scattering  agents  within  the 
cloud  does  not  vary  from  the  bottom  of  the  cloud  to  the 
top  of  the  cloud. 

These  assumptions  allow  the  two  stream  theory  and  the  diffusion  theory 
described  in  Section  3  to  be  used.  In  addition  calculated  path  lengths 
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through  oxygen  equal  the  "real"  path  lengths  in  the  reduced  atmosphere. 
Obviously,  conclusions  reached  in  the  reduced  atmosphere  must  be 
unfolded  back  to  the  real  atmosphere  before  the  task  is  completed. 


5.  2  TRANSMISSION  OF  AN  OVERCAST 

Transmission  of  a  cloud  at  a  location  associated  with  the 

area  element  at  point  p  in  the  image  plane,  Aa^,  *-s  calculated  from  the 

red  filter  picture.  The  relative  exposure  on  this  element  of  area  of 

the  film,  E  ,  is  used  in  Equation '4.  19  to  calculate  N  the  radiance  of 
P  P or 

the  cloud  element  as  viewed  at  angle  a  from  the  normal.  Given  the 
zenith  angle  of  the  sun  and  the  radiance  of  the  cloud  base, N  ,  equations 
3.  12  and  3.  14  are  combined  to  express  the  value  of  cloud  transmission. 

_  N  p  , 

T  ^ 71  a  _ ?: _  5  i 

zp  3  H©  (1+2  cos  a)  cos  z 


where 


and 


H  =  the  solar  irradiance 

9 

z  =  zenith  angle  of  the  sun 
a  =  zenith  angle  of  line  of  sight 


5.3  RATIO  OF  TRANSPORT  MEAN  FREE  PATH  TC  CLOUD  THICKNESS 


Equation  3. 12  may  be  rewritten 


m 


cos  z  +  1/  2 
T 

z 


-1 


5.  2 


This  provides  a  value  of  m  for  diffusion  theory  calculations 


5.4  OVERCAST  THICKNESS 

In  the  following  paragraphs  the  thickness  of  an  overcast,  B,  is 
calculated  from  the  base  height  of  the  overcast,  its  transmission  and  the 
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observed  transmission  of  O ^  in  the  path  from  sun,  through  overcast, 
to  detector.  Figure  5.  1  will  be  a  useful  reference  in  the  development 
of  the  working  equation.  The  equations  will  be  written  in  a  manner 
appropriate  to  the  reduced  atmosphere  (Section  3.  3  above). 

The  equivalent  sea  level  path  length  through  atmospheric  oxygen 
in  the  path  from  sun  to  cloud  top  to  cloud  bottom  to  observing  point,  Y,  is 
easily  written  for  the  reduced  atmosphere  as  Z  sec  a  +  (H-B-Z)  sec  z. 

If  r  is  the  mean  increase  in  transit  time  (in  the  reduced  atmosphere) 
for  photons  within  the  O ^  absorption  band  traveling  from  the  top  of  the 
cloud  to  the  bottom  then  the  mean  path  length  through  oxygen  from  the 
top  of  the  atmosphere  to  the  detector  is 


Y  +  c  t  =  c  t  +  Z  sec  a  +  (H-Z-B)  sec  z  5.  3 

°2  °2 

For  this  path  the  ratio  of  an  on-band  to  an  off-band  reading  using  the 
filters  described  in  Section  4.  2.  6  is 

2"  =  6  e'0-51'^"  4.27 

°  (repeated) 


or 


+  Z  sec  a  +  (H-Z)  sec  z-  B  sec  z 


5.4 


When  the  approximation  of  Equation  3.  36  holds  then  r  =  r  where  r 

U2 

is  the  mean  transit  time  of  the  cloud  at  a  wavelength  outside  the  oxygen 
absorption  band.  Now  Equation  3.  23  can  be  used  to  account  for  r  ;  that  is. 


2.  36  B 


(cos  z  +  0.  5 
T 


1.26 


) 


0.  94 


_B 

c 


(5.  5) 
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Figure  5.  1.  Cloud  in  real  (left)  and  reduced  atmosphere  (right) 


Making  this  substitution  in  Equation  5.4  and  solving  for  B  gives  the 
following  2 


B 


-  Z  sec  a  -  (H-Z)  sec  z 


Qj0.708( 


cos  z  -  0.  5 


0.  94 

-  1.  26  j  + 


5.  6 


1  |  -  sec 


The  symbol  Z  has  been  introduced  for  reasons  given  below.  In  Equation 
5.  6,  Q  3  1. 

An  effort  was  made  to  improve  on  the  approximation  t  =  r 

^2 

of  Equation  3.  36.  Many  pairs  of  curves  such  as  those  shown  in  Figure  3.  5 

were  generated  using  various  possible  values  of  Y  and  of  T.  From  the 

calculations  (which  generated  t  and  r  )  the  dependence  of  the  ratio 
—  —  °2 

of  T  to  r  Was  studied  as  a  function  of  Y  and  T.  an  empirical  relation- 
°2 

ship  was  developed  which  is 


Q  *  0.  99  x  10°-  097  l°S  T  +  IQ'0-  8  •  °-  29  (l°*  T  +  2)  log  Y 


5.  7 


This  value  of  Q  was  used  in  data  reduction.  In  retrospect  it  seems  its 
use  was  probably  not  justified  by  the  accuracy  of  the  total  data  reduction 
process. 

With  thin  clouds  and  with  significant  ground  reflection  poor  values 
of  cloud  thickness  are  obtained  because  use  of  the  value  of  T  in  equation 
5.  6  assumes  a  cloud  thick  enough  that  diffusion  theory  holds  and  also 
assumes  no  ground  reflection.  Much  of  the  data  taken  in  this  research 
effort  involved  thin  clouds  and  surfaces  which  were  the  sides  as  opposed 
to  the  bottoms  of  clouds.  To  avoid  the  errors  arising  from  these  deviations 
from  ideal  conditions,  most  of  the  data  has  been  reduced  with  cloud  thick¬ 
ness  assigned  from  best  guesses  of  the  conditions  existant  when  the  data 
pictures  were  taken. 
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5.  5  MEAN  TIME  DELAY  IN  OVERCASTS 


L 


With  a  value  of  B  calculated  using  Equations  5.  6  and  5.  7  and 
a  value  of  T  calculated  using  Equation  5.  1  the  value  of  r  can  be  deter¬ 
mined  using  Equation  5.  5.  This  yields  the  mean  time  delay  in  the 
reduced  atmosphere.  The  mean  time  delay  in  the  real  atmosphere,  r  , 


*■  T  7 


0.  708 


CH  (e-h/H-e-(h+BI/H> 


cos  z  +  0.  5 
T 


.26^ 


When  a  value  of  B  is  measured  or  assumed  then  the  following 
reasoning  is  used.  The  path  length  through  oxygen  outside  the  cloud  is 
defined  in  the  reduced  atmosphere  as  z  sec  a  +  (H-Z-B)  sec  z.  The 
total  mean  path  length  through  oxygen  is  measured  as 


In  r  -  In  6 
off 


and  the  difference  must  be  the  mean  path  length  in  the  cloud.  That  is 

r  _  .  -i  2 


In  2" .  -  in  6 
off 


Z  sec  a  -  (H-Z-B)  sec  z 


Converting  to  the  real  atmosphere  gives 


In  -  In  6 
off 


,  -<h+B)/H.  _A/1  -hi 

(e  )  sec  Z  +  ( 1  -  e 


T-  -B. 

Hc 


H,  1 

)sec  a  I 


.  -h/H  -(h+B)/  H. 

(e  -  e  1 
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5.  6  GEOMETRIC  DELAY  FROM  CLOUD  TOP  TO  SATELLITE 


In  data  reduction  a  coordinate  system  is  assigned  to  points 
on  the  data  film  using  the  edges  of  the  rectangular  frame  of  the  picture 
as  the  directions  of  the  x  and  y  axis..  The  relationships  between  coor¬ 
dinates  of  a  point  on  the  film  of  the  camera-sphere  system  and  the  azimuth 
and  zenith  angles  of  the  corresponding  line-of-sight  in  space  have  been 

discussed  in  4.  2.  1  above.  A  coordinate  system  X  ,  Y  and  Z  is  set  up 

s  s  s 

in  space  with  the  camera  atX  =  Y  =Z  =0,  with  Z  the  positive 

s  s  s  s 

vertical  and  with  X  parallel  to  x,  the  film  coordinate,  and  Y  parallel  to  y. 
s 

If  x  and  y  are  the  film  plane  coordinates  of  the  intersection  of  the 
o  o 

camera  optical  axis,  then  \  ,  the  clockwise  azimuth  of  the  point  x,  y  from 
the  positive  x  axis  is: 


for  x  positive  and  y  negative 

for  x  negative  and  y  positive  or  negative 

for  x  positive  and  y  positive 


X  =  -tan 

X  =  180°  -tan 
X  =  360°  -tan 


-1  y  '  yo 


x  -  X 


-1  y  -  yo 


o 

-1  y  -  yQ 
X  -  X 


The  clockwise  azimuth  of  the  line-of-sight  corresponding  to  point  x,y 
is  also  X.  The  zenith  angle  a  of  such  a  line-of-sight  as  shown  in  5.  1  is 


a  =  tan 


.1(^777^7 


1 


+  2 


sin  < 


D^(x-xo)2+(y-yJ2 

R  f 


— ! 


5.  11 


where  f  is  the  focal  length  of  the  camera  lens  D  is  axial  distance  from  camera 
lens  to  surface  and  R  of  sphere  is  radius  of  sphere. 

Corresponding  to  x  and  y  in  the  camera  focal  plane,  the  coor¬ 
dinates  X^,  Y^,  Zh  of  the  intersection  of  the  line-of-sight  with  the  base  of  the 


overcast  at  altitude  are 

X,  =  (h  tan  a)  cos  X 

h 

Y,  =  -h  tan  a  sin  X 
h 

Zh  -  ■> 

The  coordinates  of  the  top  of  the  overcast,  Xt>  Y  ,  Z^,  above  this  point  are 


=  Z,  +  B  =  h  +  B 
t  h 

Where  B  is  the  overcast  thickness. 

The  difference  between  the  distance  from  the  point  Xt»  Yt>  Z t 
to  the  satellite  and  the  distance  0,0,0  to  the  satellite  is  Ap  (X^.Y^.Z  ) 
given  by 

Ap  (Xt,Yt,Zt)»  -Vzt2  +  P2  cos  [a-tan-1(|-)]  g 

where  P  =  +  Y^  cos  (X  -y  +  ,3) 

7  =  azimuth  of  satellite  (clockwise  from  north) 

3  -  azimuth  of  positive  X  axis  (clockwise  from  north) 

s 
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The  difference  between  the  light  transit  time  from  X^,  Y  ,  to  the 
satellite  and  that  from  0,  0,  0  is  A  t(X^  Y  Z^). 

AP(Xt,Yt,Zt) 

At  (Xt.  Yt,  Zt)  =  - g - 


1  Jy.  2  2  -l  p. 

*  -  —  Tzt  +  Pt  cos  [a-  tan  (j-)  ]  5.13 

The  mean  time  of  flight  of  light  from  the  source  at  0,0,0  to  the  top  of  the 

cloud  at  X  Y  Z  t  (X  Y  Z  )  given  by 
t  t  t  t  t  t 

7(X.,Y  ,Z  >  =  h  sec  +  T  5.14 

where  t  is  given  by  equation  5.  10. 

The  mean  time  delay  of  light  which  starts  from  0,0,0,  leaves  the  top  of 

the  overcast  at  X  ,  Y  ,  Z  and  is  appraised  at  the  satellite  At  (a, 7)  given 
t  t  t 

by 


At  (cr,  7)  =  A  t  (X^,  Y^,  Z^)  +  h  sec  a  +  t  5.15 

5.  7  THE  RADIANCE  CENTROID 

In  the  geodetic  technique  under  study  source  location  is  determined 
by  differences  in  times  of  arrival  at  several  satellites  of  part  of  a  light  signa¬ 
ture  characteristic  of  the  source.  A  different  geodetic  technique  involves 
measurements  of  the  angles  to  the  source  from  reference  platforms  in  several 
satellites.  In  this  latter  case  the  presence  of  clouds  or  overcast  may  cause 
1  displacement  of  the  apparent  source  direction  and  an  associated  error  in 
<  4->-igned  source  location.  The  data  taken  in  the  experiments  described 
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here  can  be  used  to  specify  the  centroid  of  radiance  on  the  top  of  the 
cloud  structure  or  the  overcast.  Most  direction  sensing  equipment 
will  point  at  the  centroid. 

For  an  overcast  of  thickness  B  and  base  altitude  h  the 

altitude  of  the  centroid  Z  =  Zt  =  B  +  h.  Equations  for  the  other  two 

c  t 

coordinates,  Xc  and  Yc  are  developed  in  the  following  paragraphs. 

The  photographs  are  made  up  of  a  finite  number  of  elements 

of  equal  area  Aa  (pixels)  located  at  discrete  poisitions  x.,  y,  ,  where  j 

3  k 

runs  from  1  to  401  and  k  runs  from  1  to  401.  The  data  is  reduced  so 
that  an  exposure,  E^,  is  calculated  for  each  position  x_. , .  The  centroid 


of  exposure  xc»yc> 


is  defined  as  follows: 


51 


E x. 

.3!?— 1. 


VEjk 


and 


5.  16 


= 


LI  E  y 
j  k  jk  yk 


3 


There  is  little  physical  significance  to  this  centroid  of  exposure  on  the 
film.  It  has  been  described  because  it  is  expressed  simply. 

The  important  centroid  is  that  which  would  exist  at  the  top 
of  the  overcast  or  cloud  structure  if  the  only  light  source  were  an  isotropic 
radiator  at  the  location  of  the  camera -sphere  system.  Let  the  source  be 
a  sphere  with  a  the  area  of  its  cross  section  and  N  its  radiance  which  is 
equal  to  the  radiance  of  the  sun.  The  radiance  of  the  top  of  an  overcast  is 
written,  using  equation  3.  14 


Nlxti'V 


N  a  cos 

e  o 


.2  2 

h  sec 


08  *  (3  (1  +  2 
\14  7T  (1 

O'  ' 


cos  q>  )(1  +  2  cos  z ) 

+  B/l) 


) 


5.  17 
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where  a  is  given  by  equation  5.  11 

z  is  zenith  angle  of  sun 

h  is  base  height  of  overcast 

X  .  =  h  tan  a  cos  X 

Y_  =  -h  tan  a  sinX 
tj 


The  power,  p,  radiated  to  the  sun  from  unit  area  is 
a 


p  =  N  (X  Y  )  Q  cos  z 
a  tj  tk  a 


An  identical  expression  results  for  the  power  flowing  to  a  spherical 

receiver  at  the  camera  location  with  cross  section  a  when  the  overcast 

o 


is  irradiated  by  the  sun.  That  is 


N  ft  a  cos  z  cos  a  /o  ,,  , r 
9  9  0  (3  (1+2 

oP  aP  u2  2  V  147T  (l 

h  sec  a  ' 


1+2  cos  a)(  1+2  cos  z) ' 

Mi  +  b/l) 


Note  that  at  the  point  on  the  base  of  the  overcast  intercepted  by  the  line  of 

sight  at  zenith  angle  a  the  solid  angle  subtended  by  the  receiver  (or  source) 

located  at  the  camera  position  is  (l  ,  where 

o 


n  = 
o 


u2  2 
h  sec  a 


For  each  pixel  of  area  Aa  there  is  a  corresponding  element  of  solid  angle 
A  w  given  by  equation  4.  11  as 


Aw  = 


M2  R2 
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The  area  of  the  base  of  the  overcast  intercepted  by  this  element  of 
solid  angle  is 


,2  3 

h  sec  a 


Aw  = 


n  ,2  3 

2  h  sec 

2  2 
M  R 


If  a  ,  the  receiver  or  source  cross  section,  is  taken  to  be  the  entrance 
o 

aperture  of  the  camera -sphere  system  then  the  exposure  of  the  pixel  at 

x.  y,  is 
J  Jk 


E..  =  P..  AT  =  p  A  A  T  At 

jk  o  jk  orb  hjk 


where 


At  is  the  exposure  time 


T  is  the  transmission  of  the  camera  sphere  system 

and  the  subscript  j  and  k  are  used  to  specify  values  associated 

with  the  pixel  at  x.y.. 

1  J 

When  an  isotropic  source  of  cross  section  aQ  is  used  at  the  camera  location 

the  amount  of  power,  radiated  to  the  sun  from  an  area  AA,  ..  located  at  X,„ 

hjk  tj’ 

Y  is  by  Equation  5.  19, 


P,..  p..  AA...  =  p.,  AA,..  =  ■  -j 

tjk  m  jk  hjk  o  jk  hjk  T  At 


This  important  result  relates  the  exposure  of  the  film  at  a  particular  pixel 
to  the  amount  of  power  that  would  reach  the  sun  given  a  source  radiating  into 
the  solid  angle  corresponding  to  the  pixel.  The  source  in  this  case  has 
radiance  equal  to  the  sun  and  cross  sectional  area  equal  to  the  entrance  pupil 

i 

of  the  camera  sphere  system.  The  centroid  of  radiance  on  the  top  of  the 


overcast  is  at  X.  Y. 

tc  tc 


£  ¥  E  X 

j  k  Jk  n 

£  f  Eik 

Z  t  E  Y 

i  k  i*  Yti 
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5.  8  EFFECTS  OF  GROUND  REFLECTIVITY 


In  all  of  the  work  described  above  it  has  been  assumed  that  the 
isotropic  source  only  radiates  into  the  upper  hemisphere.  This  would  be 
the  case  if  the  isotropic  source  were  located  above  a  ground  surface  of 
zero  reflectivity.  In  the  case  that  the  ground  has  a  diffuse  reflectivity 
equal  to  p  then  reflected  light  will  contribute  to  the  power  radiated  to  space. 
If  Pracj  +  is  the  power  measured  by  all  sky  radiometer  (which  looks  into  the 
upper  hemisphere  alone) 


P 


rad  + 


g  2  2  E 

j  k 


jk 


where  g  is  a  calibration  coefficient  and  E 

x  .  y.  .  Then 
J  k 


jk 


is  the  exposure  of  the  pixel  at 


P.  .  .  =  P  .  ,  +  2g/jZ  E  cos  o 
total  rad  +  j  jk  jk 


or 


P<otal  ■  JV1  +2pcOS“'jkl 


5.  23 


The  term  cos  a  enters  because  it  is  assumed  that  the  radiation  is 
Jk 

scattered  from  the  ground  has  an  intensity  distribution  about  the  normal 
which  obeys  Lambert's  law. 
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SECTION  6 


COMPUTER  PROCESSING 

6.1  INTRODUCTION 

This  section  presents  the  computer  processing  procedures  used 
to  derive  expected  pulse  shapes  from  the  photographic  data  and  the  math¬ 
ematical  relationships  presented  in  previous  sections.  At  the  outset  a 
"hands-off"  set  of  computer  programs  was  visualized.  The  term  "hands- 
off"  as  used  here  means  a  system  in  which,  once  a  computation  is  started 
using  the  measured  input  parameters  and  the  photographic  data,  the  derived 
pulse  shape  is  a  unique,  inevitable  outcome  of  the  computations.  Such  a 
system  does  not  allow  for  human  judgement  during  the  process.  The  data 
analyst  serves  only  to  keep  the  flow  of  data  going  and  direct  the  outcome 
of  one  program  to  the  next  until  the  final  pulse -shape  plot  is  obtained. 

For  reasons  which  will  be  explained  in  more  detail  below,  a 
completely  "hands-off"  system  was  not  quite  achieved.  One  parameter, 
cloud  thickness,  proved  too  difficult  to  derive  satisfactorily  exclusively 
from  the  experimental  data  and  the  decision  was  made  to  permit  sufficient 
"hands-on"  intervention  to  supply  cloud  thickness  on  the  basis  of  typical 
clouds. 

For  either  of  these  systems  the  computer  processing  starts  with 
the  primary  photographic  data  plus  support  information  and  derives  expected 
pulse  shape.  The  final  product  of  the  system  is  a  plot  of  the  predicted, 
normalized  power  versus  time  at  the  receiver. 

Figure  6.  1  presents  the  major  features  of  the  complete  data 
processing  system  most  of  which  is  accomplished  in  the  computer.  The 
functions  named  in  several  of  the  boxes  of  Figure  6.  1  have  been  explained 
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in  other  sections  of  this  report  and  will  not  be  discussed  here.  The 
computer  processing  system  consisted  of  a  series  of  six  primary  programs 
which  were  used  sequentially,  each  providing  one  step  in  the  chain  from 
digitized  data  to  final  predicted  pulse  shape.  Several  additional  programs 
were  used  in  supporting  roles.  The  digital  data  was  organized  into  a 
square  matrix  of  data  points  centered  on  the  zenith  and  containing  the  data 
points  of  interest.  The  computer  processing  consisted  of  operating  on  each 
data  point  sequentially  and  independently  of  the  others  with  the  output  stored 
on  magnetic  tape.  All  of  the  programming  was  standard  FORTRAN  and  will 
not  be  discussed  further.  The  following  paragraphs  present  some  of  the 
problems  unique  to  the  data  being  processed  and  the  manner  in  which  these 
problems  were  handled. 

6.2  DIGITIZATION 

All  film  digitization  was  done  at  Optronics  International  Inc.  , 
Chelmsford,  MA  on  a  Model  P-1700  Photomation  Drum  Scanner.  This 
instrument  has  the  necessary  speed  to  make  a  two  dimensional  scan  of  the 
films  in  a  reasonable  time.  It  also  has  the  capability  to  write  the  digitizated 
data  on  unexposed  film  to  verify,  after  photographic  processing,  that 
digitization  has  been  properly  accomplished.  The  digitized  data  is  written 
on  9  track  magnetic  tapes  ready  for  computer  processing. 

To  prepare  for  scanning  short  strips  of  original  data  consisting 
of  four  related  exposures  were  mounted  on  a  transparent  carrier  sheet  suit¬ 
able  for  the  digitizer  drum.  Figure  6.  2  is  a  contact  print  of  six  such  strips, 
the  data  from  three  data  sets.  The  infrared  film  strips  were  all  mounted 
at  a  slight  angle  to  the  red/ blue  film  strips  to  compensate  for  a  slight 
misalignment  in  azimuth  between  the  two  cameras.  This  angle,  just  under 
2°,  was  determined  by  a  careful  examination  of  the  film  from  the  two  cameras 


in  a  projection  type  image  plane  digitizer  which  permitted  the  angle  be¬ 
tween  the  film  edge  and  lard  marks  on  the  horizon  to  be  measured.  By 
mounting  the  films  in  this  manner,  all  scans  would  follow  the  same  azimuth 
track  and  thus  avoid  image  rotation  in  the  computer  to  obtain  image  registra¬ 
tion.. 

The  scanning  parameters  chosen  were  a  50  micron  by  50  micron 
square  data  grid.  This  grid  size  is  a  compromise  between  the  spatial 
resolution  of  the  film  and  number  of  data  points  to  be  handled  per  image. 

This  grid  size  leads  to  about  12700F  data  points  in  each  circular  image. 

More  data  points  per  image  would  nave  proved  excessive  and  a  courser 
grid  can  always  be  constructed  in  the  computer  if  necessary. 

6.  3  PHOTOGRAPHIC  SENSITOMETRY 

One  of  the  two  primary  cameras  was  loaded  with  Eastman  Kodak 
Type  2476  film  for  the  red  and  blue  exposures  and  the  other  camera  with 
Eastman  Kodak  High  Speed  Infrared  film  type  HIE  for  the  ON-BAND  and 
OFF-BAND  exposures.  All  films  were  machine  processed  with  sensito- 
metric  exposure  strips  at  the  head  and  tail  ends.  Four  of  these  sensito- 
metric  strips,  a  head  and  tail  for  each  film,  type,  wore  digitized  at  the 
same  time  the  primary  data  was  digitized.  A  D-Log  E  curve  was  derived 
from  the  digitized  data,  one  for  each  film  type,  which  relates  density  on 
the  film  to  relative  exposure.  Figures  6.  3  and  6.4  are  the  final  curves  used 
in  all  subsequent  computer  processing. 

6.4  IMAGE  REGISTRATION 

It  is  an  implicit  assumption  with  this  data  processing  system  that 
the  line-of-sight  associated  with  any  given  pixel  from  one  photograph  is  the 
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same  for  all  four  photographs  for  that  data  set.  Pictures  can  fail 
to  register  because  of,  (1)  magnification,  (2)  translation  in  either  or 
both  of  two  directions,  or  (3)  rotation.  During  set-up  the  magnification 
of  the  two  cameras  had  been  adjusted  until  they  agreed  within  the  dimensions 
of  one  pixel.  Lack  of  register  in  rotation  was  handled  by  a  counter¬ 
rotation  of  the  pictures  before  digitization.  Translation  along  and  across 
the  scan  direction  is  caused  by  the  axes  of  the  two  spherical  mirrors  being 
out  of  alignment.  During  set-up  the  axes  were  aligned  by  observing 
through  the  camera  viewfinder  a  small  white  button  (the  zenith  button)  which 
had  been  accurately  glued  to  the  top  of  each  spherical  mirror  such  that  a 
ruled  cross  on  the  button  coincided  with  the  axis  of  the  sphere.  The  axes 
of  both  spherical  mirrors  were  made  vertical  with  a  sensitive  level  on  the 
mirror  baseplates. 

It  was  found  using  pictures  reconstituted  on  the  Optronics  Inter¬ 
national  P-1700  that,  when  the  zenith  buttons  of  a  picture  from  the  RED/ 
BLUE/  pair  was  superimposed  on  one  from  the  ON -BAND /OFF -BAND 
pair,  the  features  of  the  photographs  did  not  coincide  exactly.  A  trans¬ 
lation  correction  of  a  few  pixels  along  both  axes  was  required  to  bring  the 
pictures  into  register.  This  correction  was  applied  to  all  data  sets. 

6.  5  FIELD-OF-VIEW  LIMITATION 

In  Figure  6.2  it  can  be  seen  that  each  all-sky  image  occupies 
only  the  central  portion  of  the  frame.  Outside  of  the  image  area  appears 
data  board  information,  a  level,  a  watch,  and  fiducial  markers.  Prior  to 
use  as  input  data  it  is  necessary  to  transfer  to  a  secondary  data  tape  only 
those  pixels  which  pertain  to  the  all -sky  image  and  to  refer  each  image  to 
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a  common  point.  This  was  accomplished  by  locating  the  scan  number  and 
data  point  number  of  the  center  of  the  zenith  button  and  transferring  those 
data  points  which  fall  within  a  square  image  area  centered  on  the  zenith 
button  and  extending  to  a  zenith  angle  of  85°.  This  resulted  in  a  final 
grid  401  scans  wide  and  401  data  points  long,  a  total  of  160801  pixels. 

All  data  points  falling  within  a  circle  of  201  pixels  radius  were  transferred 
to  the  secondary  tape  unchanged  while  all  pixels  falling  within  the  square 
but  outside  the  circle  were  transferred  as  zero.  Figure  6.  5  is  a  typical 
image  reconstructed  from  the  data  after  such  a  transfer.  There  are  just 
under  127000  pixels  within  the  circular  image.  The  pixels  outside  the  circle 
were  retained  to  provide  scans  with  an  equal  number  of  digital  data  points, 
a  procedure  which  simplifies  the  computer  coding. 

6.  6  BLUE  SKY  FLAGGING 

The  digitized  data  from  only  one  of  the  four  primary  data  photo¬ 
graphs  is  used  to  establish  the  radiance  of  the  cloud  base,  the  data  from 
the  red  filtered  photograph.  Before  this  digital  data  can  be  use  in  a  final 
pulse  shape  calculation,  each  pixel  must  be  identified  as  one  belonging  with 
a  clear  line-of-sight  into  blue  sky  or  one  which  relates  to  a  cloud  base.  This 
was  accomplished  by  entering  the  data  from  both  the  red  and  the  blue  pictures 
into  the  computer  memory,  taking  the  blue/ red  ratio  for  each  pair  of  pixels, 
(one  blue,  one  red)  and  deciding  on  the  basis  of  the  absolute  value  of  the 
ratio  whether  the  pixel  related  to  blue  sky  or  cloud  base.  Actually  this 
was  done  a  scan  at  a  time  to  conserve  computer  memory  space.  In  a  pre¬ 
liminary  computer  run  a  histogram  was  constructed  of  the  probability  of 
occurence  for  all  values  of  blue/ red  ratio  separated  into  1000  bins.  As 


03 


expected  this  histogram  exhibited  a  doable  peak  with  very  few  of  the  ratio 
values  falling  in  the  vicinity  of  ratio  equal  to  one,  values  of  ratio  lower 
than  one  belonging  to  cloud  and  higher  than  one  to  blue  sky.  In  the  computer 
which  produced  the  input  to  the  final  calculations  the  red  radiance  was 
transferred  from  input  tape  to  output  tape  unchanged  if  the  blue/ red  ratio 
indicated  a  cloud  pixel  and  as  a  negative  radiance  if  the  ratio  indicated 
blue  sky.  In  all  subsequent  calculations  the  blue  sky  pixels  are  detected 
by  a  sign  check  on  the  radiance  value. 

As  a  check  on  this  procedure  for  blue  sky  identification  the  flagged 
digital  data  was  computer  processed  to  generate  a  two  density  mask  and  put 
the  results  on  magnetic  tape.  In  this  mask  all  blue  sky  pixels  were  given  a 
uniform  high  density  and  all  non-blue  sky  pixels  a  uniform  low  density.  When 
the  digital  data  was  written  on  photographic  film  at  Optronics  and  processed, 
the  resulting  film  becomes  an  actual  mask  showing  the  blue  sky  regions. 
Figure  6.  6  shows  the  mask  produced  by  the  blue  sky  flagging  procedure 
described  above.  This  mask  applies  to  the  data  photographed  shown  in 
Figure  6.  5. 


SECTION  7 


EXPERIMENTAL  RESULTS 

7.1  GENERAL  COMMENTS 

This  section  presents  results  obtained  from  the  experimental 
data  by  the  procedures  given  above.  The  bulk  of  the  material  in  this 
section  concerns  the  predicted  pulse  shapes  derived  from  the  cloud 
photographs.  Following  these  data  are  the  two  sections  which  present 
(1)  estimates  of  the  centroid  of  the  cloud  radiance  and  (2)  the  results 
obtained  when  an  estimate  of  the  earth's  albedo  is  included  in  the  calcu¬ 
lations. 

7.2  PREDICTED  LIGHT  PULSE  SHAPES 

The  predicted  pulse  shapes  presented  in  this  section  are 
those  obtained  by  applying  the  theory  and  procedures  given  above  to  the 
cloud  data  recorded  at  the  Kennedy  Space  Center.  All  of  the  cloud  condi¬ 
tions  analyzed  fall  into  the  general  category  of  "broken  cloud".  The 
weather  which  prevailed  during  the  six  week  period  of  data  recording 
provided  no  occurrence  of  a  solid  overcast  without  rain  (the  experimental 
equipment  was  not  intended  to  operate  in  the  rain).  This  is  indeed 
unfortunate  because  broken  cloud  is  more  difficult  to  analyze  and  less 
amenable  to  comparison  to  computer  models.  The  results  presented 
were  selected  as  the  best  representation  of  the  occasions  on  which  data 
from  all  cameras  were  available. 

The  predicted  pulse  shapes  which  follow  are  arranged  by 
date  and  time  as  the  data  was  recorded.  For  each  date  and  time  seven 
pulse  shapes  are  given  the  first  of  which  places  the  assumed  satellite 
position  along  the  line  to  the  sun's  position  when  the  data  was  recorded. 
The  remaining  six  pulse  shapes  were  derived  from  the  same  recorded 


data  but  where  the  line  of  sight  to  the  assumed  satellite  position 
has  one  of  three  elevation  angles  and  one  of  two  azimuth  angles. 

Each  group  of  seven  predicted  pulse  shapes  is  proceeded  by  a  photo¬ 
graphic  print  of  the  red  filtered  sky  data  used  to  derive  the  pulse.  The 
pulse  heights  have  all  been  normalized  to  unity. 

All  of  the  data  presented  in  this  section  with  two  exceptions 
were  derived  by  the  working  equations  which  use  an  estimated  cloud 
thickness  as  opposed  to  thicknesses  derived  from  the  data.  The  two 
exceptions  are  for  2  August  1977  at  1415  EDT  and  2  August  1977  at 
1416  EDT,  the  data  from  two  data  sets  from  the  same  data  run.  In 
these  cases  the  cloud  thickness  was  derived  from  the  data. 

Since  all  of  the  figures  in  this  section  are  self-explanatory, 
they  have  not  been  numbered  or  captioned. 


7.3  TIME  DELAYS 

The  intent  of  the  experiments  reported  here  was  to  deduce  the 
time  dependent  probability  of  delay  in  transit  of  photons  originating  from 
an  isotropic  source  on  the  ground  and  received  by  satellites  at  various 
positions  in  space.  Photons  making  direct  (unscattered)  passage  from 
source  to  receiver  would  have  zero  time  delay.  Under  these  conditions 
the  predicted  light  pulse  shapes  presented  below  would  represent  true 
delays.  Unfortunately,  the  quality  of  the  time  base  is  in  question;  however, 
errors  which  exist  in  the  absolute  time  of  delay  are  constant  for  any  set 
of  seven  curves  relating  to  one  time  of  data  acquisition.  This  allows 
relative  delays  between  lines -of-sight  to  be  compared  for  the  data  reduced 
here  because  no  direct  lines -of-sight  are  involved.  In  addition, if  the 
assumption  is  made  that  there  is  some  range  of  "reasonable"  shapes  for  the 
pulses,  then  "reasonable"  zero  times  can  be  selected  for  the  nine  sets 
of  curves  and  the  data  then  compared  to  direct  line -of-sight  timing.  This 
has  been  done.  These  concepts  are  elaborated  in  the  following  paragraphs. 
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time  was  1  second.  There  is  evidence  in  the  data  that  this  exposure  was 
sometimes  erratic  by  amounts  as  great  as  t  30  percent  of  the  average. 

The  effect  of  this  erratic  behavior  can.  at  worst,  multiply  or  divide  the 
true  on -band 'off -band  ratio  by  as  much  as  1.  7.  There  seems  to  be  no 
specific  information  in  the  data  that  can  lead  to  a  correction  of  this 
possible  error. 

One  of  the  authors  (HS)  has  inspected  each  set  of  seven  light 
pulse  shapes  and  decided  on  a  most  probable  correct  zero  time.  This  was 
done  without  reference  to  any  of  the  other  eight  sets  and,  once  selected, 
these  probable  correct  zero  times  have  not  been  changed.  These  selected 
probable  zero  times  are  believed  to  be  within  t  5  /usee  of  the  true  zero 
times.  The  selected  zero  times  are  given  in  Table  7.  1. 

Table  7.  1.  Probable  Correct  Zero  Times 


DATE 

TIME 

TIME  ON  GRAPH  OF  CORRECT  ZERO 

July  28 

1237 

-  7  /usee 

Aug  1 

1622 

-  5  /usee 

Aug  2 

1415 

+  27 

Aug  2 

1416 

+  7 

Aug  3 

0935 

+  12 

Aug  3 

0937 

+  5 

Aug  4 

1325 

-  10 

Aug  6 

1300 

-  15 

Aug  7 

1056 

0 
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The  effects  of  the  clouds  in  the  geodetic  technique  under  study 
can  be  shown  by  convolving  power  pulse  shape  characteristic  of  the  sources 
in  question  with  the  light  pulse  shapes  presented  ipro.  The  light  pulse 
shapes  shown  here  are  those  which  would  be  observed  with  an  isotropic 
delta  function  source.  Some  conclusions  of  the  cloud  effects  can  be 
reached  by  assuming  the  source  to  be  a  delta  function  and  asking  what  the 
amount  of  cloud  induced  delay  is  at  the  time  when  the  pulse  reaches  some 
fraction  of  its  maximum  value.  This  has  been  done  for  all  63  predicted 
pulse  shapes.  For  the  two  sets  of  seven  shapes  appropriate  to  3  Aug  '77 
a  value  of  60  percent  of  maximum  was  timed.  This  avoided  decisions 
arising  from  the  double  pulses  involved.  All  other  curves  were  timed  at 
80  percent  of  maximum.  The  resulting  times  are  presented  graphically 
in  Figure  7.  1.  In  this  figure  the  zero  of  time  is  that  selected  by  the 
assumption  of  reasonable  shapes  described  above.  The  crosses  (x)  repre¬ 
sent  the  delayed  time  for  the  case  of  satellite  in  the  direction  of  the  sun. 

The  altitude  and  azimuth  of  the  sun  is  given  beside  the  point.  The  circles  (o) 
represent  in  descending  order  on  the  graph  the  delay  times  for  20J,  40’,  and 
60°  elevation  at  azimuth  0°  (north)  and  60°,  40°,  and  20°  elevation  at  azimuth 
180°  (south).  The  repeatability  of  the  whole  process  carried  out  here  may 
be  judged  by  inspection  of  the  graphs  for  2  August  and  for  3  August.  For 
each  day  the  results  of  two  completely  separate  data  runs  are  presented. 

The  original  pictures  were  taken  one  minute  apart  on  2  August  and  two 
minutes  apart  on  3  August.  In  the  2  August  data  reduction  the  cloud  base 
height  was  measured  and  the  cloud  thicknesses  were  calculated  using  the 
radiance  of  the  cloud  base  as  an  input.  In  the  case  of  the  3  August  data  the 
cloud  thickness  was  assigned.  The  close  agreement  each  day  between  the 
plots  of  time  delays  for  data  taken  a  minute  or  two  apart  indicates  that  the 
overall  system  of  data  acquisition  and  reduction  is  repeatable  with  little  error. 
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7.4  CENTROID  OF  RADIAM  E 


Section  5.7  developed  the  working  equations  necessary  to 
derive  the  centroid  of  radiance  one  would  observe'  from  a  position  in 
space  if  the  cloud  conditions  which  prevailed  were  illuminated  by  an 
isotropic  source  at  the  legation  of  the  cam  >ra,  These  equations  lead 
to  the  x  and  y  coordinates  in  the  image  plane  in  pixels  (Equation  5.  16) 
or  the  X  and  Y  coordinates  on  the  top  of  the  overcast  in  kilometers 
(Equation  5.22).  Table  7.2  presents  the  results  of  the  computation  for 
six  cases. 

Table  7.2  includes  the  page  number  of  a  photograph  which 
illustrates  the  cloud  conditions  which  prevailed  for  each  of  these  results. 
The  coordinates  apply  to  these  photographs  with  upper  right  quadrant  the 
+x  and  +y  (or  the  +X  and  +Y)  quadrant  and  the  remaining  quadrants  assigned 
in  the  normal  way.  It  is  possible  to  plot  the  image  plane  coordinates  on 
the  reference  photographs  knowing  the  diameter  of  the  photographs  to  be 
401  pixels.  If  this  is  done,  great  caution  should  be  used  in  drawing  con¬ 
clusions  from  the  results.  The  photographs  illustrate  the  cloud  conditions 
that  prevailed  but,  in  no  way,  represent  the  conditions  of  cloud  illumination 
that  were  assumed  for  the  centroid  calculations.  The  photographs  were 
taken  with  the  clouds  uniformly  illuminated  by  the  sun.  The  centx-oid 
coordinates  apply  to  an  isotropic  radiation  source  located  at  ground  level 
underneath  the  clouds.  For  a  perfectly  uniform  overcast  this  mode  of 
illumination  would  result  in  the  centroid  of  radiance  appearing  directly  over 
the  isotropic  source.  It  is  not  surprising  all  of  the  centroid  locations 
given  in  Table  7.  2  fall  reasonably  close  to  the  zenith. 


Table  7.2.  Coordinates  of  Radiance  Centroid 


DATE 

TIME 
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7.  5  EFFECT  OF  EARTH  ALBEDO 

It  was  pointed  oat  in  Section  5.  8  that  the  work  presented  in 
this  report  assumed  an  isotropic  source  radiating  into  the  upper  hemi¬ 
sphere.  Equation  5.23  developed  in  that  section  was  evaluated  for  a 
single  case  to  estimate  the  magnitude  of  the  error  introduced  by  this 
assumption.  The  data  from  8/2/77  at  1414  EST  was  used  for  this  calcu¬ 
lation  using  the  worst-case  ground  reflectivity  of  unity.  The  results 
indicate  a  potential  error  of  20  percent  in  the  coupling  coefficient.  For 
more  realistic  values  of  ground  reflectivity  the  potential  error  was 
considered  negligible  in  comparison  with  other  error  sources. 
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CONCLUSIONS 


The  system  of  data  acquisition  and  reduction  developed  for 
this  project  has  produced  reasonable  results.  In  two  test  situations  data 
runs  were  made  a  minute  or  two  apart  and  were  reduced  independently. 

For  each  situation  the  associated  two  data  runs  resulted  in  almost  identical 
results.  The  data  reduction  process  was  almost  done  in  a  hands-off  manner 
however  erratic  behavior  of  the  shutter  on  the  camera  used  to  make  the 
on -band/ off -band  pictures  resulted  in  calculated  oxygen  path  length  errors 
of  as  much  as  ±  10  km.  Such  errors  were  consistent  for  all  values  derived 
for  a  particular  data  run.  A  method  for  estimating  the  magnitude  of  the 
resulting  error  was  made  based  on  the  assumption  of  a  reasonable  shape 
for  the  data -predicted  light  pulse  curves.  This  hands-on  step  is  part  of 
the  data  interpretation  process  that  was  used  in  this  report.  It  is  not 
basic  to  the  technique  developed  for  this  problem  but  it  is  used  because 
of  equipment  failure. 

This  is  a  passive  technique  in  which  steady  state  observations 
are  used  to  calculate  what  would  happen  as  a  function  of  time  to  a  light  pulse. 
The  reciprocity  theorem  was  invoked  so  that  observations  of  downward 
flowing  sunlight  could  be  used  to  predict  the  behavior  of  upward  flowing  light 
from  a  source  near  the  earth's  surface.  Since  a  pulsed  source  is  of  interest 
a  time  dependent  reciprocity  relationship  has  been  enunciated  and  has  been 
used  in  this  work.  The  validity  of  the  time  dependent  relationship  has  not 
been  demonstrated.  It  is  accepted  for  this  work  because  it  "seems 
reasonable"  and  is  required  to  get  the  job  done. 

A  single  number  has  been  chosen  to  represent  the  effects  of 
clouds  on  the  transfer  of  a  pulse  of  radiation  from  a  source  near  the  ground 
to  a  satellite.  The  number  is  the  time  in  microseconds  necessary  for  the 
delayed  pulse  to  build  up  to  some  selected  fraction  of  its  maximum  value. 
Values  for  the  fraction  of  maximum  used  here  are  0.  8  and  0.  6.  The  times 
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so  determined  are  taken  to  be  representative  of  the  cloud  induced  errors 
in  the  geodetic  timing.  For  the  63  calculations  made  in  this  report  the 
average  of  the  times  of  delay  of  the  0.  8  (or  0.  6)  of  maximum  point  is  about 
15  jusec.  There  are  seven  delay  values  for  each  of  the  nine  data  runs.  The 
maximum  spread  of  the  delay  values  has  been  determined  for  each  data 
run.  The  average  of  these  nine  spread  values  is  10  ^usec.  Figure  7.  1 
gives  a  detailed  presentation  of  ths  information. 

In  Figure  7.  1  the  plots  for  1  August  shows  two  points  with 
negative  delay.  This  is  unreasonable  and  could  be  corrected  by  a  new 
choice  of  the  zero  correction.  When  the  probable  corrected  zero  times 
of  Table  7.  1  were  established  it  was  decided  that  there  would  be  no  second 
guessing.  The  1  August  plots  result  from  this  effort  to  stay  honest. 

An  expected  correlation  between  the  half  width  of  the  time 
delay  curve  and  the  coupling  coefficient  was  not  demonstrated.  This  is 
shown  in  Table  8.  1.  The  coupling  coefficient,  CC,  is  the  ratio  of  the 
radiation  from  the  whole  sky  incident  on  an  isotropic  detector  to  the 
radiation  that  would  have  been  incident  from  the  sun  outside  the  atmosphere. 

The  radiance  of  the  base  of  overcasts  or  clouds  is  dependent 
in  part  on  the  reflectivity  of  the  ground.  This  effect  has  not  been  included 
in  the  data  reduction.  It  can  be  of  considerable  importance  especially 
when  there  is  foliage  on  the  ground.  This  is  due  to  the  rapid  increase  in 
reflectivity  of  foliage  with  wavelength  increase  between  0.1  fj  and  0.  8^. 
Future  observations  should  include  careful  study  of  ground  reflectivity. 

An  ideal  location  for  measurements  involving  absorption  at  0.  7  6,u 
would  be  at  sea. 

The  most  uncertain  measured  input  to  the  data  reduction  was 
the  height  of  the  base  of  the  overcast  or  cloud  layer.  For  data  reduction 
at  any  particular  time  a  single  value  of  base  height  was  used  for  all  clouds 
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Table  8.  1  Comparison  of  the  coupling  coefficient  with  characteristics  of 
relative  amplitude  vs  delay  curves. 


DATE 

TIME 

TIME  OF  MAX 

WIDTH  AT  1/2  MAX 

CC 

7/28 

1237 

6  jusec 

30 

0.  99 

8/1 

1622 

10 

25 

0.  57 

8/2 

1415 

38 

32 

0.  22 

8/2 

1416 

20 

27 

0.  22 

8/3 

0935 

20 

35 

43 

1.02 

8/3 

0937 

12 

43 

0.98 

8/4 

1325 

8 

30 

0.  58 

8/6 

1300 

10 

20 

0.  70 

8/7 

1056 

10 

30 

42 

0.  64 
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and  overcast  present.  This  value  was  the  average  of  several  stereo 
determinations . 


Thickness  of  cloud  or  overcast  was  calculated  using 
Equation  5.  6  for  the  data  taken  on  2  August  at  1415  and  1416.  Equation 
5.  6  assumes  that  only  the  cloud  base  is  viewed  and  it  is  horizontal.  In 
addition  it  is  assumed  that  the  cloud  is  thick  enough  to  be  described  by 
diffusion  theory.  As  a  result  of  this  sunlit  and  thin  portions  of  the  over¬ 
cast  had  extraordinary  thicknesses  assigned  to  them.  A  major  effort 
was  invested  in  attempts  to  improve  the  calculation  of  cloud  thickness 
but  in  the  end  it  was  decided  to  assign  cloud  thickness  from  any  reasonable 
input  available  (cloud  type  etc.  ).  Tests  showed  that  the  relative  amplitude 
vs  delay  time  curves  are  reasonably  insensitive  to  the  assigned  cloud 
thickness. 

The  cime  delay  curves  are  generated  from  radiance  values 
appropriate  to  the  clouds  and  overcast.  The  blue  sky  is  neglected  inten¬ 
tionally.  A  more  mature  data  analysis  technique  should  include  contribu¬ 
tions  from  the  clear  sky. 

The  principles  of  the  technique  would  not  be  modified  if  the 
data  were  taken  properly  by  some  electro-optical  system.  If  this  were 
done  there  could  be  great  saving  in  data  manipulation  before  data 
reduction  and  an  almost  real  time  system  might  be  possible. 
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APPENDIX  A 


EXPERIMENTAL  EQUIPMENT 

The  experimental  phase  of  this  project  spanned  a  six  week 
period  during  July  and  August  of  1977.  The  equipment  was  designed 
and  constructed  at  HSS  Inc  in  Bedford,  MA.  and  transported  to  Kennedy 
Space  Center,  FL.  for  installation.  The  site  provided  was  a  NASA 
photographic  station  not  currently  in  use.  The  coordinates  of  this  station 
have  been  given  in  Section  5  of  this  report.  This  site  provided  a  concrete 
pad  of  sufficient  size  flanked  on  three  sides  by  a  paved  area.  Most  impor¬ 
tant  this  site  provided  an  almost  clear  line-of-sight  in  all  directions  down 
to  a  zenith  angle  of  85  .  There  were  a  small  number  of  local  structures 
that  did  block  an  insignificant  portion  of  the  total  field-of-view  above  the 
85  3  horizon.  This  site  was  the  primary  site  and  was  designated  as  the 
A  site.  Approximately  1  km  from  the  A  site  a  secondary  site,  B,  was 
located  for  stereo  picture  purposes.  The  B  site  was  visible  from  the  A 
site  and  the  connecting  line  was  used  as  a  reference  line  for  orienting  the 
camera  formats.  The  A  site  and  B  site  were  connected  by  field  telephone. 
The  instrument  complement  at  each  site  was: 

A  Site 

Red/ blue  all  sky  camera 
On/ off  band  all  sky  camera 
Flight  Research  stereo  camera 
Automax  stereo  camera 
All-sky  radiometer 
Teleradiometer 
Sun  photometer 
Camera  controller 
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B  Site 


Flight  Research  stereo  camera 
Automax  stereo  camera 
Camera  controller 

Figure  A1  is  a  general  view  of  the  A  site.  The  two  pipe  stands  in 
the  foreground  support  the  primary  data  Nikon  cameras  and  their  reflecting 
spheric  mirrors,  the  camera  on  a  side  arm  above  and  the  mirror  on  a  side 
arm  below.  Just  in  back  of  these  cameras  are  the  A  site  stereo  cameras, 
a  Flight  Research  Multidata  4C  camera  with  60  3  field -of-view  lens  on  the 
left  and  an  Automax  camera  with  Fairchild  All-sky  lens  on  the  right. 

Both  cameras  are  mounted  in  watertight  enclosures  and  adjusted  to  bring 
the  optic  axis  of  the  lens  vertical.  The  60  0  field-of-view  of  the  Flight 
Research  camera  was  accommodated  by  a  flat  window  in  the  top  of  the 
camera  enclosure.  The  all-sky  lens  of  the  Automax  required  a  hyper- 
hemispheric  dome  to  accommodate  the  full  field-of-view.  The  all-sky 
radiometer  appears  immediately  behind  the  right  pipe  stand.  The  tele- 
radiometer  is  located  in  a  tripod  at  the  right  of  the  picture. 

Figure  A2  shows  the  two  stereo  cameras  located  at  the  B  site.  The 
cables  lead  to  the  camera  controller  located  in  a  watertight  box  off  to  the 
side.  The  square  box  sitting  near  the  base  of  the  Automax  camera  enclo¬ 
sure  is  simply  an  aluminum  foil  covered  cardboard  box  placed  over  the 
top  of  the  camera  enclosures  when  pictures  are  not  being  taken  to  reflect 
sunlight  and  keep  the  enclosure  as  cool  as  possible  for  film  protection. 

Figure  A3  is  a  detailed  view  of  the  reflecting  spheric  mirror  used  with 
on-off  band  primary  camera.  The  red/ blue  primary  camera  had  an 
identical  mirror.  The  mirror  mount  provides  for  translation  in  two 
directions,  rotation  about  the  pipe  stand,  and  accurate  leveling.  These 
motions  permit  the  accurate  location  of  the  center  of  the  small  white 
button  on  the  top  of  the  mirror,  the  zenith  button,  directly  under  the  optic 
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axis  of  the  camera  located  on  the  side  arm  above.  The  camera  itself 
was  mounted  in  a  manner  which  permitted  the  optic  axis  of  the  lens  to 
be  rotated  to  the  vertical. 

Once  the  camera  support  stand  was  set  up  and  lagged  to  the 
concrete  pad  the  primary  camera  alignment  procedure  was  as  follows. 

First  the  camera/mirror  separation  was  adjusted  to  the  desired  dis¬ 
tance.  Then  the  camera  axis  was  made  vertical  using  a  sensitive  level 
(0.  17  degrees  tilt  per  0.  1  inch  of  bubble  movement)  mounted  on  the 
camera.  The  mirror  position  was  then  adjusted  until  the  zenith  button 
appeared  centered  on  the  crosshairs  of  the  camera  view  finder.  When 
in  position  a  series  of  test  exposures  were  made  and  processed  The 
text  exposures  from  both  the  on/off  band  and  the  red/ blue  cameras  were 
measured  to  determine  the  image  distance  for  the  two  fiducial  marks 
appearing  on  the  mirror  support  plate  (see  Figure  A3).  (During  manu¬ 
facture  the  distance  between  the  fiducial  marks  on  each  mirror  support 
plate  were  made  equal.)  The  measured  distance  between  marks  for  each 
camera  was  converted  to  the  magnification  for  that  camera  and  the  two 
compared.  The  camera/ mirror  distance  for  one  camera  was  then  changed 
and  the  alignment  repeated  until  the  camera  magnifications  agreed  within 
a  resolution  element  on  the  film. 

Figure  A4  shows  a  detailed  view  of  the  all  sky  radiometer.  The 
reflecting  spheric  mirror  can  be  seen  in  the  center  of  dome.  A  collar 
surrounds  the  mirror  as  a  field -of- view  limiter.  The  reflected  energy 
enters  the  optical  system  located  in  the  tube  directly  above  the  mirror 
and  falls  on  a  silicon  detector.  The  detector  output  signal  is  carried  oy 
cable  to  the  electronics  located  in  die  housing  b'low  the  dome.  This  housing 
has  a  removable,  watertight  transparent  plate  (removed  in  Figure  A4) 
which  allows  the  comp’ete  unit  to  operate  unattended.  In  this  mode  the 
output  signal  is  recorded  on  u  strip  chart  i  "corder.  The  arm  which  appears 


in  the  photograph  is  used  to  block  the  direct  rays  from  the  sun  under 
some  conditions  of  measurement. 

Figure  A5  shows  the  teleradiometer.  This  instrument  is 
very  similar  to  the  all  sky  radiometer  except  that  it  has  very  limited  field 
of  view.  The  radiometer  housing  was  designed  to  accept  any  of  the  filters 
used  in  the  primary  data  cameras,  red,  blue,  on  band,  or  off  band.  The 
output  signal  is  read  on  the  all  sky  radiometer  electronics.  The  tripod 
mount  for  the  teleradiometer  also  supported  a  boresighted  elbow  telescope 
for  identifying  the  point  being  viewed  by  the  teleradiometer. 

The  sun  photometer  (not  shown)  is  a  hand-held  instrument 
designed  by  Fred  Volz.  The  instrument  measures  the  irradiance  of  the  sun 
in  a  narrow  field-of-view.  Once  calibrated  the  reading  obtained  can  be 
related  to  the  transmission  of  the  atmosphere  in  the  spectral  band  being 
observed  from  observer  to  outside  the  atmosphere.  Each  instrument 
provides  several  wavelength  bands  of  measurement. 

The  experimental  data  recording  procedure  was  as  follows  using 
a  three -man  crew,  two  at  A  site  and  one  at  B  site.  When  a  cloud  structure 
of  interest  appeared,  the  B  site  operator  was  notified  by  telephone  to  turn 
on  the  camera  control  units,  make  the  stereo  camera  identification 
exposures,  and  stand  by  for  a  count.  The  same  procedure  was  carried 
out  at  the  A  site.  The  identification  exposures  were  made  on  both  stereo 
cameras  and  served  to  identify  the  next  group  of  stereo  pictures  as  to  site, time, 
date,  camera,  and  film  number.  When  both  sites  were  ready  a  verbal 
count  was  started  over  the  phone  to  synchronize  the  stereo  pairs.  At 
the  A  site  the  two  operators  immediately  made  the  four  primary  data 
exposures.  B  site  was  then  again  called,  the  count  repeated  and  more 
stereo  pairs  exposed.  The  four  primary  data  exposures  were  then 
repeated  followed  by  a  final  group  of  synchronized  stereo  pairs.  This 
sequence  completed  a  data  recording  set.  The  operation  of  the  all  sky 
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radiometer,  teleradiometer,  and  sun  photometer  were  not  included 
in  a  primary  data  set  run.  These  instruments  were  operated  in  a 
support  mode  to  collect  calibration  and  additional  data. 
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